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Abstract: We consider the infrared structure of hadron-hadron collisions at next-to-next- 
to leading order using the antenna subtraction method. The general form of the subtraction 
terms is presented for double real, real-virtual and double virtual contributions. At NLO 
and NNLO it is shown that the virtual and double virtual subtraction terms can be written 
in terms of integrated dipoles, formed by systematically combining the mass factorisation 
fvj contributions and integrated antenna functions. The integrated dipoles describing I unre- 

solved partons, denoted J 2 , are related to Catani's IR singularity operators, l\- {e). It 
is shown that the IR pole structure of the virtual and double virtual contributions can be 
written as a sum over integrated dipoles within the antenna subtraction formalism and 
the master expressions analogous to Catani's one- and two-loop factorisation formulae are 
derived. To demonstrate the techniques described in this paper, we apply antenna subtrac- 
tion to the production of two gluon jets via quark- antiquark scattering at NLO and NNLO. 
Double real, real-virtual and double virtual subtraction terms are explicitly derived for the 
leading colour NNLO contribution. 
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1. Introduction 



In QCD, the (renormalised and mass factorised) inclusive cross section for a hard scatter- 
ing process in proton-proton collisions that produces a particular final state, X, has the 
factorised form, 

where the probability of finding a parton of type i in the proton, carrying a momentum 
fraction £, is described by the parton distribution function [£p)d£ and the partonic 
cross section d&ij for parton z to scatter off parton j and produce X, normalised to the 
hadron-hadron flux 1 is summed over the possible parton types i and j. As usual and 
Hf are the renormalisation and factorisation scales which are frequently set to be equal for 
simplicity, fiR = fxp = A 4 - 

For suitably high centre of mass scattering energies, the infrared- finite partonic cross 
section has the perturbative expansion 



a s {Vn)\ A ~NLO , ( a s{HR)\ 2 A ~NNLO , 



where the next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) strong 
corrections are identified. In general, the leading-order cross section will be proportional to 
an overall power of the strong coupling, a™, with the value of m depending on the process 
under consideration. NLO computations for multiparticle processes are becoming standard 
for comparison with LHC data. However, for processes with only one or two objects in 
the final state, the spectacular performance of the LHC and LHC experiments means that 
NNLO accuracy is mandatory. 

Calculations of higher order QCD corrections to the underlying Born-level process 
encounter additional integrals that must be combined before physical predictions can be 
extracted: loop integrals in the case of virtual corrections and phase-space integrals for 
real corrections, reflecting the additional particles present in higher order calculations. It 
is well known that for theories containing massless states, both the virtual and real radiative 
corrections are peppered with IR singularities which conspire to mutually cancel to form 
the finite physical cross section. 

In order to achieve this cancellation, the IR singularities must be extracted from the 
real emission phase space integrals in the form of a Laurent expansion in the dimensional 
regularisation parameter e, where D = 4— 2e. Subtraction schemes are a well established so- 
lution to this problem and several methods for systematically constructing general subtrac- 
tion terms have been proposed in the literature at NLO [1-5] and at NNLO [6-24]. So far, 
successful applications of subtraction at NNLO to specific observables were accomplished 
with antenna subtraction [9], (^-subtraction [25], sector decomposition [26-29], non-linear 
transformations [19], and most recently with an approach based on sector- improved residue 



lr rhe partonic cross section normalised to the parton-parton flux is obtained by absorbing the inverse 
factors of £i and £2 into d&ij . 
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subtraction [18,20]. The c/T-subtraction method is restricted to processes with colourless 
final states at leading order and is based on the universal infrared structure of the real 
emissions, which can be inferred from transverse momentum resummation. This method 
has been applied at NNLO to the Higgs production [30], vector boson production [31,32], 
associated VH- production [33], photon pair production [34] and in modified form to top 
quark decay [35]. The sector decomposition approach relies on an iterated decompostion 
of the final state phase space and matrix element, allowing an expansion in distributions, 
followed by a numerical evaluation of the sector integrals. It has been successfully ap- 
plied to Higgs production [21, 36, 37] and vector boson production [38] at NNLO. The 
sector-improved residue subtraction extends NLO residue subtraction [2] , combined with a 
numerical evaluation of the integrated subtraction terms. It has been applied to compute 
the NNLO corrections to top quark pair production [39,40]. 

The antenna subtraction formalism that is the focus of this paper has been fully de- 
veloped for all types of scattering processes at NLO [9, 41-43] and further developed at 
NNLO [9,44-47]. Within this method, the subtraction terms are constructed from so-called 
antenna functions which describe all unresolved partonic radiation (soft and collinear) be- 
tween a hard pair of radiator partons. The hard radiators may be in the initial- or final- 
state, and in the most general case, final-final (FF), initial-final (IF) and initial-initial (II) 
antennae need to be considered. The subtraction terms and therefore the antennae also 
need to be integrated over the unresolved phase space, which is different in the three con- 
figurations. Recently, all of the integrals relevant for processes at NNLO with massless 
quarks have been analytically evaluated [44-47]. 

The antenna subtraction approach has been successfully applied to the infrared struc- 
ture of three-jet events at NNLO [48-50] and the subsequent numerical calculation of the 
NNLO corrections to event shape distributions [51-55], the moments of event shapes [54,55] 
and jet rates [53,56]. For hadron colliders, the antenna method has been applied at NNLO 
to the all gluon contribution to dijet production [57-59] and the production of heavy par- 
ticles [60-63]. 

There is however another aspect to the use of the antenna subtraction method. The 
antenna approach successfully isolates the infrared singularities, which themselves have a 
very particular structure. Therefore, one can ask the questions: 

1. How does the infrared structure guide the construction of the real radiation subtrac- 
tion terms? 

2. How do the integrated subtraction terms relate to the known infrared structure of 
the loop amplitudes? 

The main goal of this paper is to develop the systematics of the antenna subtraction scheme 
at NLO and particularly at NNLO by focussing on the structure of the subtraction terms 
for the double real, real-virtual and double virtual channels. 

It is well known that the explicit poles of the virtual contributions are described by 
Catani's one- and two-loop factorisation formulae [64, 65] and follow a colour-dipole struc- 
ture, reflected by the IR singularity operators 1^' '(e) and which describe the IR 
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singularities associated with two colour connected hard particles. It is similarly well estab- 
lished (and in fact the keystone of all subtraction methods) that real radiation amplitudes 
factorise in IR divergent limits following an antenna factorisation pattern whereby a pair of 
hard partons radiate unresolved partons. The unintegrated antenna functions are used to 
mimic the implicit divergence of the real contributions while the poles associated with the 
integrated antenna functions directly cancel the explicit poles of the virtual contributions. 

The implicit divergence of the NLO and NNLO cross sections is captured by distinct 
blocks of subtraction terms which follow a predictive structure based on the colour ordering 
of the matrix elements under consideration. After integration and combination with the 
relevant mass factorisation contributions, a number of structures emerge which we call 
integrated dipoles. At one-loop, J^p is related to an integrated three-particle tree-level 
antenna and describes the unresolved radiation between two colour connected particles. 

(2) 

Similarly, at two- loops, J 2 fulfills the same role for double unresolved radiation and 
involves integrated four-particle tree-level antennae, three-particle one-loop antennae and 
products of three-parton tree-level antennae. 

The IR structure of any two-loop contribution is given by one- and two-loop integrated 
antenna strings which are simply formed from and J^' '. For a given colour ordering 
of a particular n-particle process, j\% ^ and j\% ^ can immediately be written down as a sum 
over integrated dipoles that involve colour connected particles, 



where £ = 1,2 denotes the single and double unresolved integrated antenna strings re- 
spectively and the sum is over colour connected pairs of partons. Equivalently, the pole 
structure is identified as a sum of dipole-like terms proportional to (\sij\) ~ le that link the 
colour connected particles i and j. 

Of course, the form of J$ and jff in terms of integrated antennae imposes a partic- 
ular structure on the unintegrated antennae that make up the real radiation subtraction 
terms. Understanding the explicit pole structure of virtual amplitudes in terms of inte- 
grated antenna strings gives a direct connection between the block structure of the unin- 
tegrated subtraction terms and the explicit pole structure of virtual contributions, thereby 
simplifying the construction of the double real, real-virtual and double virtual subtraction 
terms. 

The paper is organised in the following way. In Section 2 the general structure of the 
antenna subtraction terms relevant for NLO calculations is discussed, in both unintegrated 
and integrated forms. This procedure leads to a natural definition of the NLO integrated 
dipole Section 3 then extends this analysis to NNLO where we construct groups of 

antenna subtraction terms that remove the double and single unresolved singularities from 

(2) 

the double real and real-virtual contributions. The NNLO integrated dipole J 2 naturally 
emerges from the groups of integrated antenna functions and mass factorisation kernels 
and (along with combinations of J^) reproduces and properly subtracts the explicit poles 
of the double virtual contribution to the NNLO cross section. 




(1.2) 
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As an example of the power of the structures motivated in general terms in Sections 
2 and 3, Section 4 considers the production of gluon jets from quark-antiquark scattering 
at NLO and NNLO. We derive the double real, real- virtual and double virtual subtraction 
terms relevant for the NNLO corrections to the qq — > gg channel, which is one of many- 
channels that contribute to the dijet rate at the LHC. In Section 5 we revisit the double 
virtual subtraction term for e + e~ — > 3 jets at NNLO and show that it has the same 
structure in terms of integrated antenna strings. Finally, our findings are summarised in 
Section 6. For completeness a number of appendices are enclosed detailing the collinear 
splitting kernels in the N, Np basis we employ and giving an explicit expression for the 
integrated soft factor with an initial-final mapping. 



2. Antenna subtraction at NLO 

The structure of NLO calculations within the antenna subtraction formalism is particu- 
larly simple. Nevertheless it demonstrates some of the key features relevant to NNLO 
calculations which are fully explained in this section. 

We are interested in the higher-order corrections to the underlying leading-order (LO) 
cross section, which for the hadro-production of n jets is given by, 

daijMZiHi'b^) = [ d«7? LO (a#i,6#2), (2.1) 

J n 

where the Born-level cross section is obtained by evaluating the tree-level contributions 
for (n + 2)-parton scattering with partons i, j in the initial-state carrying a fraction £i ; 2 
of the parent hadron's momentum R\^- This quantity is then integrated over the final- 
state phase space where we keep track of the number n of final state particles using the 
shorthand notation, 



/ 

J n 



(2.2) 



Explicitly, the Born-level cross section is given by, 

&Pij,LO = NlO Yl d$ n(P3,-- - ,Pn+2-,Pl,P2) ^~ 
X M° +2 (<7(1,--- ,71 + 2)) 4 n) (Mn) + 



1 

AT2 



(2.3) 



where S n is a final-state symmetry factor, and {p} n denotes the set of (n) final-state 
momenta. M® +2 denotes an (n + 2)-parton squared partial amplitude for a given colour 
ordering denoted by a and i labels the parton with momentum p^. The factor Mlo w dl 
in general contain all non-QCD factors and some overall QCD factors such as the overall 
power of the coupling. The 2 — > n particle phase space is defined by, 



d$„(p 3 ,--- ,p n+2 ;pi,P2) = 

\d x d 



d j-i p3 d d - l Pn+2 
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The jet algorithm builds n jets from m final-state partons with momenta labelled by 
the set {p} m - In the case of the leading-order cross section defined in Eq. (2.3) m = n and 
there is an exact parton-jet correspondence. 

The NLO correction to the n-jet cross section contains three contributions and is given 

by, 

d&ij,NLO = / d°ij,NLO + [ d °Yj,NLO + d °ij,NLo)i ( 2 - 5 ) 
Jn+1 Jn\ J 

where d&fj and AaJ- NLO are the real and virtual NLO corrections and dafj^LO ^ s ^ n e NLO 
mass factorisation contribution, 

d^Lo(6tfi,6#2) = "/ (°^) C(e)rM M ( Zl ,z 2 ) dauMzi^z&H,), 

(2.6) 

where (7(e) = (47r) e e~ e7 . The mass factorisation contribution serves to remove all initial- 
state collinear singularities from the cross section by absorbing them into the redefined 
physical Parton Distribution Function (PDF) and can be written in terms of the one-loop 
Altarelli-Parisi kernels, 

r SiU*i. *2) = *(1 - z 2 )S lj T^(z 1 ) + 6(1 - z 1 )5 kt T < i f(z 2 ). (2.7) 

The real cross section is known to contain soft and collinear IR divergences and so a 
subtraction term, dfrfj jvlO' ^ s constructed from antenna functions and reduced multiplicity 
matrix elements. This subtraction term must remove all implicit singularities from the real 
emission cross section without introducing further spurious singularities of its own and 
be local in the sense that the subtraction is successful point-by-point in phase space. A 
secondary requirement of a subtraction term is that it be analytically integrable in order to 
convert the implicit divergence of the subtraction term into the explicit singularities of the 
integrated subtraction term. 2 When combined with the mass factorisation contribution, 
these terms conspire to cancel the explicit poles of the virtual contribution and produce a 
finite NLO cross section. The requirement of analytic integrability is automatically satisfied 
in the antenna subtraction method as all three- and four-parton antenna functions have 
been successfully integrated for FF, IF and II configurations. 

Given these considerations, the NLO cross section can be reorganised in the following 
way such that each square bracket is free from both implicit divergence and explicit poles, 

daij,NLO = [d°ij,NLO ~ daf j>NLO ] + / [da(^ NLO - daf j>NLO ] , (2.8) 

Jn+l Jn 

where the virtual subtraction term is given by the real subtraction term integrated over 
the single unresolved phase space and the NLO mass factorisation contribution, 

dvfj^LO = ~ daf j NLO - daf-Jj L0 . (2.9) 



2 One could equally well numerically integrate the subtraction term as in Ref. [18] or [23,24]. 
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2.1 Construction of the real emission subtraction term 



The single real emission cross section takes the form, 

M?j,NLO = N$LO d $n+l(P3,--- , Pn+3', Pi , P2) 



M n ° +3 K1, ■ ■ ■ ,n + 3)) 4 n+1) (Mn + i) + ^(^) 



n+1 

(2.10) 



where, as in Eq. (2.3), S n+ \ is a final-state symmetry factor, {p} n+ \ denotes the set of 
(n + 1) final-state momenta, M^ +3 denotes an (n + 3)-parton squared partial amplitude 
for a given colour ordering denoted by a and i labels the parton with momentum pi. The 
overall coupling is given by 



^^o(^)gg (2.11) 



where C(e) = C(e)/8vr 2 



2.1.1 Colour ordered strings 

Before discussing the various subtraction terms, it is useful to review the possible structures 
available for colour ordered matrix elements. There are two basic objects from which all 
colour structures can be assembled: gluon strings and quark strings. A gluon string is 
characterised by a trace over SU(N) generators, where au denotes the generator associated 
with external gluon k, e.g., an n-gluon colour ordered matrix element is associated with 
the colour structure, 

Xr ( T ai T a 2 ... T a n _ 1T a^ _ ( fll a 2 • • • a n _ia„) . (2.12) 

The colour ordered partial amplitude associated with this colour structure contains IR 
divergences between colour connected partons, i.e., those external legs whose associated 
generators are adjacent in the colour structure. Due to the trace nature of the colour factors 
the partial amplitudes display cyclic symmetry, which means that the colour connection 
not only exists between adjacent partons in the notation defined above, but also between 
the endpoints a± and a n . 

The other basic colour structure is the quark string, given by a string of gluon gen- 
erators bookended by fundamental quark indices, e.g., a quark string containing a quark- 
antiquark pair and n gluons has the colour structure, 

( T «i T «2 . . . T a n - 1T a n) _ (oifl2 . . . an __ ian)ij} (2.13) 

where i,j denote the fundamental (antifundamental) colour indices carried by the quark 
(antiquark). The IR divergences of the partial amplitude associated with this colour struc- 
ture also exist for unresolved configurations involving colour connected partons. The gluons 
at each end of the quark string are not colour connected to each other but are colour con- 
nected to the quark endpoints instead. 
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If the quark-antiquark pair become collinear to form a composite gluon, the colour 
structure pinches down to a trace structure and the surviving reduced matrix element 
(multiplied by the universal splitting function) tends to a gluonic partial amplitude, 

(aia 2 ■ ■ ■ a n -ia n )ij ^ {a x a 2 ■ ■ ■ a n _ia n a (ij -)). (2.14) 

Colour structures involving multiple quark-antiquark pairs are built from multiple 
quark and gluon strings where the various colour ordered amplitudes are distinguished 
by the ordering of the gluons within a quark string. Sub-leading colour structures also 
arise at the amplitude level according to the distribution of fundamental indices among the 
endpoints. For example, the four-quark two-gluon partial amplitudes have twelve colour 
structures. Six of the structures are given according to the permutations of gluons within 
the structures, 

(aia 2 )u(0)kj {ai)u{a 2 )kj {0)u{aia 2 ) kj , (2.15) 

where i, k and j, I denote the fundamental and anti-fundamental indices of the quarks and 
antiquarks respectively and {0)u = 6u denotes an empty quark string whose colour factor 
is simply a Kronecker delta in fundamental indices. The six sub-leading colour structures 
are obtained by the substitution jol, 

As discussed above, in the case of a quark string, a collinear limit between the quarks 
causes the quark string to collapse onto a gluon string. With multiple quark strings, 
a quark-antiquark collinear limit can cause quark strings to merge, with the composite 
gluon contributing to the gluon string by connecting the fundamental indices of the quark 
endpoints, e.g., 

Ki • • • a i n )ij( a ji ■ ■ ■ a j m )ki Q ^ Ki ■ ■ ■ a-i n a(kj)a jl ■ ■ ■ a jm ) u . (2.16) 

The behaviour of the colour structures in various unresolved limits is mirrored by the 
resulting reduced partial amplitudes. When multiple strings are present in a single colour 
ordered partial amplitude, the partons involved in each quark are separated from those 
involved in other strings by a semi-colon, e.g., 

M n{--- ; Qi, 9i, 92, qi;q3, 9i, &,■■■)■ (2.17) 
2.1.2 Sub-leading colour 

The sub- leading colour contribution cannot, in general, be written in the form of an in- 
coherent sum of squared partial amplitudes. However for low multiplicity final-states the 
sub-leading colour contribution can often be rewritten as the square of a coherent sum of 
QCD partial amplitudes where one or more gluons behave in an Abelian fashion. 

Matrix elements containing "Abelian gluons" obey the usual QCD factorisation for- 
mulae in all the unresolved limits of the non-Abelian gluons. The Abelian gluons do not 
couple to the non-Abelian gluons and only couple to quarks; therefore they can only be 
considered colour-connected to the quarks. In this way a sub-leading colour matrix ele- 
ment can be considered to contain two colour structures: a pure QCD colour structure 
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formed from all quarks and non-Abelian gluons, the factorisation properties of which are 
identical to those of leading colour squared matrix element, and a QED-like colour struc- 
ture containing quark pairs and colour disconnected Abelian gluons [48]. For example, a 
squared partial amplitude with two quarks, n non-Abelian gluons and one Abelian gluon, 
i, is denoted M® +3 (q, g±, • • • ,gi, - ■ ■ ,g n ,<j) and schematically has the colour structure, 

(q, 1, ■ ■ ■ , v ■ ■ , n, q)ij ~ (q, 1, ■ ■ ■ , n, q)ij <g> (q, i, q) {j . (2.18) 

The Abelian gluon can have collinear limits with the quarks and also become soft, 

M%+3fa 9u ■ ■ ■ ,9i,--- ,9n,q) — P qg ^ Q (z) M° +2 (Q, 5 i,--- ,flf„,g), (2.19) 

M n+zi.^g^--- ,9n,q) ^ S qi q M% +2 (q , g 1 , ■ ■ ■ ,g n ,q), (2.20) 

where P qg ^Q(z) and are the time-like quark-gluon splitting function and soft function 
respectively. The collinear limit with the antiquark is given by exchanging q -H- q. When the 
sub-leading colour matrix elements can be re- written in terms of squared matrix elements 
involving Abelian gluons then the IR divergent limits can be subtracted using antenna 
functions just as at leading colour. For six or more coloured particles, the sub-leading colour 
contribution cannot be written purely in terms of squared matrix elements with Abelian 
gluons. However, because the IR divergent limits of the remaining contribution are known, 
antenna subtraction can be readily applied. Notwithstanding these complications, the sub- 
leading colour contributions will in general contain fewer and less intricate divergent limits 
than their leading colour counterparts. Any finite contribution to the cross section requires 
no subtraction and is integrated numerically. The remainder of this section will focus on 
the leading colour contribution for clarity but it is clear that a similar treatment can be 
reproduced for sub-leading colour contributions. 

2.1.3 Subtraction terms 

The NLO subtraction term is based on the single unresolved limits of the relevant matrix 
elements. In our case, the subtraction term reflects the antenna factorisation of the squared 
partial amplitude. 

The full subtraction term is a sum of contributions of the type, 



d°NLO = A/jvLO 52 52 d$ n+ i(p 3 ,--- ,p n +3-,Pl,P2) ~7T~ 

— iJn+1 
perms j 

x X$(;j,.) M„° +2 (... ,j,-..) 4 n H{p}n). (2.21) 



where X® is the NLO antenna function, M® +2 the reduced matrix element with a particular 
colour ordering and the sum runs over all the possible unresolved particles j. There are 
three separate cases that correspond to the hard radiators being in the final state (FF), 
both in the initial state (II) or one in the initial-state and one in the final-state (IF). For 
DIS processes daff^Q = while for e + e~ annihilation da^ 1 ^ = da^Lo = ® m ^ NLO it 
is possible to write the subtraction terms explicitly for all configurations at leading colour 
and for an arbitrary number of partons. 
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When the unresolved parton j is colour connected to the final-state hard radiators i 
and k, the subtraction term for the partial amplitude M° +3 (- • • ,i, j, k, • • • ) takes the form, 

M° +3 (- •• ./•..,./.-.•••) • X 3 °(*, j, k) M° +2 (- ■ ■ , 7, K, ■ ■ ■ ). (2.22) 

The species of antenna function is determined by the identities of the partons contained 
within the antenna, as shown in Tab. 1. The sum over j takes into account all the possible 
unresolved partons in the colour-ordered matrix element fitting the final-final configura- 
tion. The final-final phase space map [66], (i,j,k) — > (I,K), ensures that the momenta 
involved in the antenna function are mapped onto two hard composite momenta. The IR 
divergence associated with the configuration where j becomes unresolved is described by 
the appropriate antenna function and, in the singular limit, the subtraction term tends 
to the value of the real emission cross section. The various antennae and reduced matrix 
elements appropriate for a particular colour ordered matrix element in the real emission 
are listed in Tab. 1. There are five possibilities reflecting the different particle assignments 
and possible colour structures. 



Final-Final Unintegrated Antennae 


Matrix element, M° , 3 


Antenna, X 3 


Reduced matrix 
element, M^ +2 


(" " " ! iqijgi kq\ • • • ) 


A° 3 (i,j, k) 


(••• ;Iq,Kq,---) 


( ' ' ' ! iq , jg , kg , ■ ■ ■ ) 


d° 3 (i,j,k) 


(••• ;Iq,K g ,---) 


(' ' ' j iq'i jq, kq, ■ " " ) 


Ei(i,j,k) 


(• • • ; Iq',K g , ■■■) 


(' ' ' ' igiJgi kgi ■ • • ) 


m,3,k) 


(••• ,Ig,K g ,---) 


(' ' ' ' igijq'i kq, • • • ) 


G° 3 (i,j,k) 


(••• ,Ig,K g ,---) 



Table 1: The NLO antennae, X®, and reduced matrix elements, M° +2 , appropriate for the various 
particle assignments and colour structures in the real radiation partial amplitudes, M° +3 for the 
final-final configuration. 

In configurations with initial-state partons, an initial-final antenna is necessary and the 
initial- final phase space map [43] — > (1,7) is employed to generate the composite 
momenta for the reduced matrix element so that 

M° +3 (--- X,i,3, ■■■) ^X 3 a ^(i,i,j) M n ° +2 (- «/,•••) (2.23) 

where the initial-state parton labelled l a where a is the species of the initial-state parton, 
either q or g for a quark or gluon respectively. The various antennae and reduced matrix 
elements appropriate for a particular colour ordered matrix element in the real emission 
are listed in Tab. 2. In hadron-hadron collisions, an analogous subtraction term exists for 
the second initial-state parton, i.e., l a — > % where b is the species of the second initial-state 
parton. For processes where a = b, we drop the subscript a — > b for simplicity and only 
retain the label in the species changing cases. 

The last configuration that is allowed for hadron-hadron collisions is the initial-initial 
configuration where an unresolved parton is emitted between two initial-state partons of 
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Initial-Final Unintegrated Antennae 


Matrix element, M° +3 


Antenna, A3 


Reduced matrix 
element, M® +2 


(" " " j lq> igi jqi ' ' ' ) 




(' ' ' ! Igj Jqi ' ' ' ) 


{'''■> !<?' igi 3gi " ') 


4iUij) 


(" " " 1 Iff' ' " " ) 


(" " " j iqijgi Iff' ' " " ) 


d° 3 (i,j,i) 


(' ' ' 1 Jqi Iff' ' ' ' ) 


(' ' ' i Iff'' iq'i jqi ' ' ' ) 


E° 3 (li,j) 


(" " " j !<?') " " " ) 


{■ ■ ■ i Iff, igi jgi - - - ) 


fi(U,j) 


(• • • , Iff, ' ' ' ) 


(" " " i^-giiq'i jqi ' ' ' ) 


G° 3 (U,j) 


(• • • , Iff, Jgi ' ' ' ) 


V ' iiqi^-gi jqi ' ' ' ) 


-<4,g^ q (i,hj) 


(' ' ' i Iff' «^ff' ' ' ' ) 


(• • • ; 2<j> igi Jgi ■ ■ • ) 


-ds^qii'hj) 


(' ' ' i^-qi Jgi ' ' ') 


(' ' ' i ^q' , lq! Jq, " " " ) 


U) 


(' ' ' 1 Jq' 1 Iff 5 ' ' ' ) 


(" " " iigi^-q'i jqi ' ' ' ) 




(• - - ) </j, Iff) - - - ) 



Table 2: The NLO antennae, A3, and reduced matrix elements, M° +2 , appropriate for the various 
particle assignments and colour structures in the real radiation partial amplitudes, M^ +3 for the 
initial-final configuration. The identity changing antennae are collected at the bottom of this table. 

species a and b. The subtraction term uses an initial-initial antenna function and the 
appropriate initial-initial phase space map [43] (1, i, 2) — > (1,2), 

M n ° +3 (. • • , 1 , i, 2 b , ■ ■ ■ ) — ► A 3 ° a ^^ d (i, i, 2) M n ° +2 (- --XX,--- )• (2-24) 

The various antennae and reduced matrix elements appropriate for a particular colour 
ordered matrix element in the real emission are listed in Tab. 3. As in the initial-final case, 
the subscripts are only retained in the species changing cases, a ^ c or b 7^ d. 



Initial-Initial Unintegrated Antennae 


Matrix element, M° +3 


Antenna, X 3 


Reduced matrix 
element, M^ +2 


(' ' ' i !<?' igi 2<?> " " " ) 


4j(U2) 


(■ ■ ■ j 1<?) 2^; • • • ) 


(' ' ' i Iff' igi 2ff> • ' ' ) 


£)§(i,»,2) 


(' ' ' ' Iff' 2g, • • • ) 


(' ' ' > -^ff' igi 2ff) ' ' ' ) 


f$(U§) 


(' ' ' ' -^ff' 2ff, ■ ■ ■ ) 


(' ' ' i 1?' 2<j, iq', • • • ) 


— ^3,ff-5>g(l' 2> «) 


(■ ' ' 1 2^; • • • ) 


(' ' ' ' *9' Iff' 2ff) " " " ) 


-4^,1; 2) 


(" " " ' 1?) 2g, ■ ■ ■ ) 


(" " " ) ' 2gj iqi " ' " ) 


2, i) 


(' ' ' > Ig'j 2 ff , ■ ■ ■ ) 


(■ ■ ■ 'Iff' 2g~; ^<j, • • • ) 




(' ' ' ' Iff' 2ff) " " " ) 



Table 3: The NLO antennae, A3, and reduced matrix elements, M° +2 , appropriate for the various 
particle assignments and colour structures in the real radiation partial amplitudes, M° +3 for the 
initial-initial configuration. 
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2.2 NLO mass factorisation term 

The notation for the NLO mass factorisation contribution was presented at the beginning of 
this section in Eq. (2.6). Writing the Born cross section out in terms of the matrix elements 
being integrated over the final-state phase space the mass factorisation contribution at NLO 
is given by, 

Kr v sr^ / \ f dzi d^2 1 

d(T ij,NLO = -N NLO 2^ d ®n{P3,-- ■ ,Pn+i;Xip 1 ,X 2 P2) J 



perms 



x rgjU*!, * 2 ) M°(- • • ,k, ■ ■ ■ , l, ■ ■ • ) jW(p 3 , • • • , Pn+ i). (2.25) 

In this formula there is an implicit sum over k and I, denoting the particle types of the 
initial-state partons. The explicit forms for the various mass factorisation kernels and their 
associated colour decompositions are listed in Appendix A.l. The overall factor is given 
by, 

rv _ \ r (ol s N\ R 



Klo = ^lo ( C(e) = A/jvlo C(e). (2.26) 

2.3 Construction of the virtual subtraction term 

The virtual subtraction term is constructed from the integrated real subtraction term 
and the mass factorisation contribution. For a single colour connected pair in the colour 
ordering of the virtual matrix element, a virtual subtraction term can be constructed of 
the type, 

d °NLO = -MnLO I —— -i- V d$ n (p 3 ,--- ,Pn,X 1 p 1 ,X 2 P2) 
J X 1 X 2 tin permg 

x 4\l,K) M„° +2 (... ,I,K,-..) 4 n) ({p}n). (2.27) 

I,K 



As illustrated in Fig. 1, the singular integrated dipole factor, is directly related to 
the integration of the A3 antenna function over the single unresolved phase space, and is 
defined to all orders in e. The precise form of depends on the species of the resolved 
hard radiator partons in the real radiation matrix elements that produced this integrated 
dipole contribution and the kinematic configuration, i.e., final-final, initial-final or initial- 
initial. The full virtual subtraction term is formed by joining together the various integrated 
dipoles to form a single integrated antenna string for each reduced matrix element. 

In the final-final configuration, the generic virtual subtraction term for a pair of colour 
connected partons /, K is produced by integrating the subtraction term associated with 
the real radiation matrix element M® +3 (- • • ,i,j,k,---) with the mapping (i, j, k) — > (I, K) 
and has the form, 

M° +3 (... ,i,j,k,---) -^4\l,K) M° +2 (... ,/,#,■■■)■ (2.28) 

This equation is the analogue of Eq. (2.22). The details of the integrated dipoles depends 
on the species of parton under consideration and all possibilities are listed in Tab. 4. 
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i I 




Figure 1: NLO virtual structure. Integration of the tree-level real emission over the single unre- 
solved phase space for particle j generates an integrated dipole (I, K) for each dipole pair /, K, 
containing the explicit IR poles of the virtual contribution. 



Final-Final Integrated Antennae 


Matrix element, M® +3 


Integrated dipole, J% 


Reduced matrix 
element, M° +2 


(' ' ' ! jgi kq', • " " ) 


4 1 \l q ,Kg)=A° 3 (siK) 


(' ' ' j Iqi Kq'i ' ' ' ) 


( ' ' ' 'i^qi jgjkg, ■ ■ ■ ) 


J i i\l q ,K g ) = \Vl{s IK ) 


(" ' ' i Iqi Kgt ' ' ' ) 


(' ' ' j iq'i jq'i kqi ' ' ' ) 




(' ' ' ! Iqi Kgi ' ' ' ) 


( ' ' j ig-i jgi kgt ' ' ° ) 


4 1 \l g ,K g ) = 1 ^(s IK ) 


(■■■ ,Ig,K g ,---) 


(' ' " i ^gijq'i kqi ' ' ' ) 


J% F {h,Kg) = Gl{s IK ) 


(' ' ' , Ig, K g , • • • ) 



Table 4: The correspondence between the real radiation matrix elements, M^ +3 and the integrated 
NLO dipoles and reduced matrix elements, M° +2 for various particle assignments and colour 
structures for the final-final configuration. 

When one of the partons making up the hard dipole is in the initial-state then the 
analogous formula to Eq. (2.28) now contains the appropriate mass factorisation kernels 
from Eq. (2.25), 

M° +3 (- • • , i, ■ ■ ■ ) — > J«(l, J) Af° +2 (- • • , 1, J, • • ■), (2.29) 

where the only modification is the form of the integrated dipoles, listed in Tab. 5. In this 
equation the subtraction term involving parton 1 is presented; in hadron-hadron collisions 
the initial-final virtual subtraction term involving the other initial-state parton is obtained 
by the substitution 1 — > 2 

The final kinematic configuration is when both hard radiators in the dipole are in the 
initial-state. In this case the relevant subtraction term is, 

M n ° +3 (- , 1,^,2,---) — > ^ 1} (1,2) M° +2 (... ,1,2,...), (2.30) 

where once again, only the form of the integrated dipole is different from Eqs. (2.28) 
and (2.29) and the initial-initial integrated dipoles are listed in Tab. 6. 

2.4 Summary 

The relationship betweem the real radiation subtraction term da s and the virtual sub- 
traction term da T is shown in Fig. 2. The full virtual subtraction term is generated by 
summing over the colour connected pairs in the virtual matrix element's colour ordering. 
For example, consider an (n + 3)-parton squared partial amplitude, containing a quark- 
antiquark pair and (n + 1) gluons, that maps onto the (n + 2)-parton squared partial 
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Initial-Final Integrated Antennae 


Matrix element, M^ +3 


integrated dipole, J 2 


Reduced matrix 
element, M° +2 


(" " " j igi jqi ' ' ' ) 


4 1) (i„j (7 ) = ^ f/ ( SiJ )-C(x 1 )5 2 


(' ' ' j Igi "^g! ' ' ' ) 


(■ ■ ' ; lg, igi jgi " " " ) 


4 1) (i„j fl ) = ipO (7 ( SIJ )-rW(x 1 )«5 2 


(" " " j 1?) Jgi " " " ) 


(••• ;* g ,ifl,i fl> ---) 


4 1} (jg, = ^3%, w (sij) - MVofc 


(■■• ;Jg,l g ,---) 


(" " " j l<j' i iq'i jqi ' ' ' ) 


J 2,N F \\^g) = 2^3,<7',gg( s lj) 


(" ■ ■ j Igj Jgi ' " " ) 


(' ' ' ' Iff) J's> ' ' ' ) 


__ 1 n / \ 1 -m— 1 ( 1 ) / \ 1- 

J^(l fl , J fl ) = - \ v gg^2 


(' ' ' > Jgi " ') 


(" " " i Igi *gi Jg> " " " ) 


J 2,N F 0-gi J g) = ^3,<y( s lj) 


(" " " j lg; ^ffj " " " ) 


(" " " 1 *g> Jg! " " " ) 


4!^g(lg ; J«) = -\A,gA S u) ~ 4ffVl)<*2 


(' ' ' i lg> ' ' ' ) 


(' ' ' ! *g> Igijgi ' ' ' ) 




(" " ' j Igj ' " " ) 


(' ' ' j > lg! Jg> " " " ) 




(' ' ' ! Jqi Igj ' ' ' ) 


(" " " j 2g J lg> Jgj " " " ) 


4!^g(^ 5 %) = -QUsu) - r#(*i)fc 


(' ' ' ' ^ff> Iff) ' " " ) 



Table 5: The correspondence between the real radiation matrix elements, M° +3 and the integrated 
NLO dipoles and reduced matrix elements, Af° +2 for various particle assignments and colour 
structures for the initial-final configuration. For brevity 6(1 — Xi) = Si for i = 1,2. 



Initial-Initial Integrated Antennae 


Matrix element, 


Integrated dipole, 


Reduced matrix 
element, M° +2 


(" " " j lg> igi *2q, • • • ) 


j^(t q , %) = Al Q ,(s u ) - r%( Xl )h - C(x 2 )5 1 


(' ' ' ) lg> 2g! • • • ) 


(" " " ! 1?) igi 2 S , • • • ) 


J$\l q , %) = Vl Qq (s u ) - F q 1 q ) (x 1 )S 2 - \T a 1 a \x 2 )5 1 


(' ' ' ) lg> 2ff) " " " ) 


(• • • j Iff) ^ff) 2 ff , ■ ■ ■ ) 


J 2 1] (ig,%) = J$ M (* H ) ~ irff (xijft - £r#(x 2 )*i 


(' ' ' ) Iff) 2ff> • • • ) 


(" " " ) 1?) 2g,iq] • • • ) 


4!U(lg)2g-) = -^,gff( S 12) -Tgff (^1 


(' ' ' ) lg> 2g! " " " ) 


(■ ■ ■ j ^g) Iff) 2 ff , ■ ■ ■ ) 


4!U(lg)2ff) = -^,gg( s i2) " rgff J (xi)5 2 


(" " " j 1?) 2<j) ■ ■ ■ ) 


(' ' ' ) lg' ) 2gJ 4g, " " " ) 


■^(Mff) = -^ 3 °a'a,«( s i2) " r£Wi 


(' ' ' ! lg') 2ff) " " " ) 


(' ' ' > Iff) 2^-; ig, • • • ) 


4!g^(lff)2ff) = -^(a H ) - rJJWi 


(' ' ' ) Iff) 2g 5 " " " ) 



Table 6: The correspondence between the real radiation matrix elements, M" +3 and the integrated 
NLO dipoles and reduced matrix elements, M° +2 for various particle assignments and colour 
structures for the initial-initial configuration. For brevity 5(1 — x\) = 5\, 5{1 — X2) = 62- 



amplitude M® +2 (l q ,i a , • • • ,j g ,2 q ). Once the leading colour single unresolved subtraction 
term for this matrix element has been integrated and combined with the mass factorisation 
terms, the resulting virtual subtraction term is simply given by, 



^nlo = -Nnlo j ^1^1 1- V* d$ n (p 3 ,--- ,p n ;x 1 p 1 ,x 2 p2) 
J x\ x 2 b n ^—^ 



perms 



x JnUhiig,--- ,jg,%) M Q n+2 (\ q ,i g ,--- ,j g ,2 q ) jW(p 3 , • • • ,p n ), (2.31) 
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da T : 


da T ~ J^MO 












> 






da s : 


d«7 S ~ X 3 M„ 



Figure 2: The NLO correction to n-parton scattering requires the construction of the subtraction 
term dafj LO to remove all IR divergent behaviour of the (n + l)-parton real emission. Integrating 
this subtraction term and combining with the NLO mass factorisation terms generates the virtual 
subtraction term, d<r^ L0 . Antenna subtraction allows da]^ LO to be generated by first constructing 
d^NLO or vice versa - 

where the full integrated antenna string is formed from the sum of integrated dipoles, 

^n+2(^5' igiJgi 1 kg, lg, 2q) = J 2 ^ ig) + J 2 ^ i}gi 3g) + ' ' ' 

+ J i 2 1 \k g ,l g ) + J2 1 \%,l g ). (2.32) 

Each term in J^\ 2 ls proportional to a particular kinematic factor, (\sij\)~ e . We there- 
fore see that the singularities present in each term correspond to a particular singular 
contribution iff- The full singularity structure is simply obtained by summing the inte- 
grated dipoles to form or equivalently by summing Using Catani's one-loop 
factorisation formula [64], 

Voles{Ml +2 {l,--- ,n + 2)) = 21%(e; 1, • • • , n + 2)M° +2 (1, • • • ,n + 2), (2.33) 

where I^\ 2 { e 'i 1> ' ' ' > n + 2) is the n-parton IR singularity operator, which can be written 
as a sum of dipole contributions, 

/« 2 ( e ;l,---,n) = ^jSj ) (e; % ), (2.34) 

and the sum runs over colour connected pairs of partons. It is clear that the virtual sub- 
traction term, formulated in terms of integrated antenna strings, reproduces the explicit 
poles of the virtual cross section to complete the IR finite calculation at NLO. It is noted 
that now the correspondence between the poles of ano - ^ ne one -l° P matrix elements 
has been established, the explicit pole cancellation can be carried out without recourse 
to 1^+2 • This correspondence is particularly clear because of the way the integrated an- 
tenna string has been constructed as a sum of dipole-like antenna contributions and by 
recognising that the mass factorisation contributions naturally fit into the structure of the 
integrated antenna strings so that only the genuine final-state IR poles remain in the vir- 
tual subtraction term. This puts the virtual subtraction term in a particularly convenient 
form for carrying out the explicit pole cancellation against the virtual contribution. 
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3. Antenna subtraction at NNLO 



The NNLO cross section involves two additional particles in the scattering process, ap- 
pearing as either real or virtual particles. These additional particles lead to three distinct 
contributions: double real da^^ NLO , real- virtual da^ v NNLO or double virtual do - ^ jvjvlo> 
containing (n + 4)-, (n + 3)-, (n + 2)-partons and zero-, one-, two-loops respectively, where 
n is the number of jets in the final-state, 

&Oij,NNLO = / &°fj R NNLO+ I [^fj^NNLO + ^ijNNLo) 
Jn+2 Jn+1 v 7 

+ / {^YjfNNLO + d ^,NNLo) > t 3 - 1 ) 
Jn 

where J d$ is defined in (2.2). The NNLO mass factorisation terms are naturally parti- 
tioned into terms containing (n + 1)- or n-parton final-states, which are given by, 

J J Z\ Z 2 

x (°^) C(e) r% kl ( Zl , z 2 ) (da^ NLO - da s khNL0 ^j (z&H^ z 2 &H 2 ), (3.2) 



dz\ dz2 
zi z 2 



2 



C(ef A(z 1 ,z 2 )da kl , LO (z 1 ^H 1 , z 2 ^ 2 H 2 ) 



a s N 

+ (^) C(e) r^ kl ( Zl ,z 2 ) (dal hNLO - da T khNL0 ^ (z^H,, ^ 2 ff 2 )}, (3.3) 

(2) 

where the new ingredient at NNLO, r^ -/ fcJ , can be written in terms of one- and two-loop 
Altarelli-Parisi kernels, 

F^Uzuz,) = 5(l-za)^T2>(zi) + «(^^^ (3.4) 

Using the antenna subtraction method, subtraction terms can be constructed from 
antenna functions (now including four-parton antennae to describe double unresolved ra- 
diation and one-loop three-parton antennae to describe the single unresolved radiation of 
one-loop matrix elements) to remove all implicit divergence of the double real and real- 
virtual contributions. Introducing a double real subtraction term dcr?- jvtvlo' a real-virtual 
subtraction term daYj S NNLO , their corresponding integrated forms and the NNLO mass fac- 
torisation contributions, the NNLO cross section can be reorganised such that each term 
in square brackets is rendered free from implicit divergence and explicit poles and thus 
suitable for a numerical implementation. 

d(?ij,NNLO = / [dvfjjfNLO ~ daf j NNLO ] 

Jn+2 



+ / [do-fj V NNLO — dafj NNLO ] 

Jn+1 
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+ / [d&ijJfNLO ~ ^ij,NNLo] > 



(3.5) 



n 



where the real-virtual and double virtual subtraction terms are given by, 




(3.6) 



(3.7) 



Note that because integration of the subtraction term da^NLO & yes contributions to both 
the (n + 1)- and n-parton final states, we have explicitly decomposed the integrated double 
real subtraction term into a piece that is integrated over the single unresolved particle 
phase space and a piece that is integrated over the phase space of two unresolved particles, 



The construction of subtraction terms at the double real, real-virtual and double vir- 
tual levels requires an intimate understanding of the explicit IR singularity structure and 
implicit IR divergent behaviour of both the physical matrix elements and the subtraction 
terms for different final-state multiplicities. The structures of these subtraction terms are 
dictated by two main considerations: the colour connection of partons and explicit IR pole 
cancellation. The former completely informs the construction of the double real subtrac- 
tion term as there are no explicit IR poles to consider; the latter is the sole concern of 
the double virtual subtraction term which is free from implicit IR divergences. The real- 
virtual subtraction term contains both explicit poles and implicit divergent behaviour and 
so its structure must take into account both parton colour connection and explicit pole 
cancellation. 

As we will discuss in Sees. 3.1, 3.2 and 3.3 below, the dafj NLO , daJj NLO and da^f NLO 
subtraction terms each contain several components that each fulfill a particular function 
in constructing the correct infrared structure for NNLO computations. It is the nature of 
any subtraction scheme, that whatever is subtracted from a final state with a particular 
multiplicity, must be added back in either integrated or unintegrated form to a final state 
with a different multiplicity. To illustrate the connections between the various terms, 
Fig. 3 shows how the subtraction terms da^ NLO , daJj NLO and da^ NLO are linked. The 
"roadmap" shows how the various terms are related. The five terms introduced in the 
double real subtraction term are directly linked with specific terms in the real-virtual and 
double virtual subtraction terms. Similarly, terms in the real-virtual subtraction term are 
composed partly of integrated terms from the double real subtraction term j l daf- ^at^o 
and new virtual subtraction terms, daJ- s NNL0 . The remaining terms, f 2 dafj nnlo ano - 
Ii d&Yj S NNLO ' combine to produce the structures present in dafj NNLO ■ As the discussion 
unfolds, we anticipate that it will help the reader to refer back to Fig. 3. 

3.1 Construction of the double real subtraction term 

At NNLO the double real radiative (n + 4)-parton correction to the pp —> n jets process, 
must be included integrated over the full (n + 2)-parton phase space, subject to the con- 




(3.8) 
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straint that n jets are observed. In order for this to be achieved, the IR divergent behaviour 
of the partonic cross section must be isolated. In the antenna subtraction formalism this is 
achieved by constructing a subtraction term from three- and four-parton antenna functions 
and reduced multiplicity matrix elements. 

By construction the double real subtraction term mimics the IR divergence and fac- 
torisation of the real radiation matrix elements in all relevant single and double unresolved 
limits. At NNLO there are four colour configurations to consider when constructing the 
subtraction term for n-jet production [9,48,50,57]: 

1. d<T 5 ' a 

A single unresolved parton but the remaining (n + 1) final-state partons form n jets. 

2. da s > b 

Two colour-connected unresolved partons, i.e., two unresolved partons radiated be- 
tween a single pair of hard radiators. 

3. d<r 5 ' c 

Two almost colour-connected unresolved partons, i.e., two colour disconnected unre- 
solved partons sharing a common radiator, and including large angle soft radiation. 

4. da s > d 

Two colour disconnected unresolved partons, i.e., two colour disconnected unresolved 
partons with no radiators in common. 

Here we have suppressed the dependence of the subtraction terms on the initial-state 
partons i and j. The four possible colour configurations for the single and double unresolved 
limits of the (n + 4)-parton matrix elements provide a natural way to divide the subtraction 
term into pieces, as has been emphasised in previous works [9,48]. It should be noted that 
although the factorisation of the matrix elements is strictly classified according to the colour 
connection of unresolved partons, the associated subtraction terms readily communicate 
with one another in most unresolved limits due to the existence and cross cancellation of 
spurious singularities. These systematic cross cancellations suggest that structures exist 
in the subtraction term other than those dictated by the colour connection of unresolved 
partons. 

In the following sections we make a further reorganisation of the double real subtraction 
term into five contributions, 

d°NNLO = ^NNLO + ^NNLO + ^NNLO + ^NNLO + ^NNLO' ( 3 - 9 ) 

This manifestation of the subtraction term is broadly organised along the lines of colour 
connection but is also informed by considering how the double real subtraction terms appear 
in real-virtual and double virtual subtraction terms after integration. Such an organisation 
of the double real subtraction term allows a transparent understanding of how various terms 
from the double real level cascade down the calculation upon integration. Understanding 
this structure also permits a systematic construction of the double real subtraction term 
in terms of predictable blocks of terms. 
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3.1.1 Single unresolved subtraction term, d&ff^ LO 

The removal of single unresolved limits from the (n + 4)-parton contribution to the re- 
jet cross section at NNLO strongly resembles the subtraction term constructed to isolate 
single unresolved limits from the (n + 4)-parton contribution to the (n + l)-jet cross section 
at NLO. The two subtraction terms differ only in the number of jets allowed by the jet 
function, 



dd NNLO (Mn+4) = da^ LO ({p} n+4 ) 



An+l) j(n+l) 
J n+1 ~^ J ri 



(3.10) 



The detailed discussion about constructing the NLO subtraction terms given in Sec. 2 
is then sufficient to construct the single unresolved NNLO subtraction term d&fj^ L0 . 
Although this part of the subtraction term removes single unresolved divergences of the 
physical matrix elements for n-jet selecting observables, it also generates additional spurious 
singularities in the almost colour connected and colour disconnected limits, over-subtracting 
the divergence in each case. This is because the (n + 3)-parton reduced matrix elements 
allow an additional parton to become unresolved and yet still form n-jets. This contribution 
is reintroduced at the real-virtual level upon analytic integration where they cancel the 
explicit poles in the real- virtual contribution. They form a distinct set of terms proportional 
to the (n + 3)-parton matrix elements. 

3.1.2 Four-parton antenna subtraction term, d&ff£ LO 

At NNLO there are essential new ingredients in the form of four-parton tree level antenna 
functions that are required to faithfully reproduce the colour connected double unresolved 
divergences of the physical matrix elements. The momentum map associated with a four- 
parton antenna function is the (n + 4) — > (ra + 2) map that maps the four antenna momenta 
i,j,k,l down to two composite momenta I,L, which is different for the final-final [66], 
initial-final and initial-initial [43] configurations. This group of terms has the form, 

d&NNLo = Nnnlo Y] d* n+2 (p3,...,j> n +4;pi,P2) -J— 

x Y, X° 4 (i,j,k,l) M n ° +2 (... ,I,L,...) 4")({p}n), (3-11) 

where N^^LO contains the overall QCD coupling, colour and non-QCD factors appropriate 
for the double real radiation contribution, 

* f NNLO=X LO [-^) cgp. (3-12) 

This term is generic to all three kinematic configurations, (FF, IF and II), with the specific 
antenna functions and momentum maps depending on which configuration the term belongs 
to. The four-parton antenna functions display many different types of divergent behaviour 
including single unresolved and almost colour connected singularities all of which must be 
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properly removed elsewhere in the subtraction term. The analytic integration of the four- 
parton antenna functions is carried out over the double unresolved antenna phase space 
and so the terms in daff^ LO are reintroduced as part of the double virtual subtraction 
term. 

S b 

3.1.3 Four-parton single unresolved subtraction term, da^^j^Q 

As mentioned in the previous section, four-parton antenna functions contain single un- 
resolved spurious singularities which must be removed to ensure a proper subtraction of 
the IR divergence in the physical matrix elements. For each four-parton antenna the sin- 
gle unresolved limits are removed by constructing a subtraction term along the lines of 
d&NNLO now a PPly m S this method to the four-parton antenna function rather than 
the physical matrix elements. As such, the terms are built from three-parton antennae 
(used to remove the single unresolved limits) multiplied by another three-parton antenna 
(the remnant of the four-parton antenna after the single unresolved limit is taken) and a 
reduced matrix element which in the single unresolved limit maps on to the matrix element 
associated with the four-parton antenna. For a four-parton antenna function this block has 
the generic form 



d& NNLO = -NnNLO J2 d<Z> n+2 (,P3, ■ ■ ■ ,Pn+4-,Pl,P2) 

x Y, X°(iJ,k) X$(I,K,l) M° +2 (... ,/',£,•■■) 4 n) (Mn), (3-13) 



where the sum is over all partons in the four-parton antenna which admit a single unresolved 
singularity. As with other parts of the subtraction term, this block of terms may contain 
spurious almost colour connected singularities of its own which can arise when the secondary 
antenna in (3.13) also contains single unresolved limits. The block of terms associated with 
each four-parton antenna encapsulates all of the single unresolved singularities of the four- 
parton antenna. As discussed in Sec. 3.2, a collection of terms which correctly mimic the 
single unresolved limits of the four-parton antenna functions performs a specific role when 
reintroduced as part of the real-virtual subtraction term. 

3.1.4 Almost colour connected subtraction term, der^y^Q 

The almost colour connected contribution only exists when there are at least five partons 
in the scattering process. It is intimately linked to the four-parton antenna functions that 
contain almost colour connected unresolved partons, denoted by X®; X® antennae do not 
give a contribution to da^ c NL0 . The classification of the four-parton antenna functions 
into X® and X® types is displayed in Tab. 7. 

S b 

In almost colour connected limits, the terms in da^^ LO tend to over-subtract the diver- 
gences of the associated X®. In the same limits da^ a NLO contributes twice the subtraction 
required by the matrix elements. Both of these over-subtractions have to be accounted for 
by da^f NLO , which also includes the wide angle soft subtraction term [48,67]. 

The terms appropriate for du^ c NLO can be generated in the following fashion. We 
consider a final-final double unresolved configuration for clarity but the strategy also applies 
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Configuration 


xl 


^4° 


Final-final 


A®, B®, C4 , D 4jQ , i?4 j0 , £4, ^4 i(l , G4, G4, H® 


^4' u k,ci ^4,6' r 4,fe' 


Initial-final 


ao rO /oo ^<o r<o ttO 

t - / 4 ! Lt 4) ^4) -"4 


^4 ;-°4 -> E i >-^4 


Initial-initial 


40 R o (~<o n° f° r 10 r° n° 

Jl 4 , tf A , O4, ^ 4iad j^4,adj' ^4' -°4 


^4' -^4,n.adj ' ^4 ' ^4,n.adj 



Table 7: The classification of the four-parton antenna functions into those containing almost 
colour connected limits or not. The final state D®, E\ and F° antennae are decomposed into 
sub-antennae [48, 57] for numerical implementation. 



j 1 




Figure 4: The three regions associated with the radiation of the primary unresolved parton. Region 
I: between the hard radiators of the four-parton antenna, j, I. Region II: one place to the left in 
the colour ordering between partons a,j. Region III: one place to the right in the colour ordering 
between partons I, b. 

to initial- final and initial- initial configurations. Each X®, in daffy LO generates a block 
of terms in da^ c NLQ . For example, consider the antenna X2(j,i,k,l) with partons i, k 
unresolved and j, I the hard radiators, such that the reduced matrix element is given by 
M°(. . . , a, (jik), (ikl),b, . . .). To construct daff NLO we first consider the colour ordering 
with i and k removed, i.e., (• • • ,a,j,l,b,---) as depicted in Fig. 4. We then allow partons 
i and k to be radiated from this underlying ordering; i is radiated first, followed by k (plus 
the reverse ordering). The second radiation (k) always takes place between the antenna's 
hard radiators (region I) whilst the first radiation (i) may be inserted directly between 
(region I) and one place either side (regions II and III) of the antenna's radiators in the 
colour ordering, with a relative minus sign for these two contributions. Written in terms 
of antennae this contribution has the form, 

+ l - Xi(j, i, I) X$((ji), k, (if)) MX . . , a, (55fc), (£$0), b, . . .) 

- \ X° 3 (a, i, j) Xl((ji), k, I) M°(. . . , (oT), ((ji)k), (kf),b,.. .) 

- \ Xi(l, i, b) Xi(j, k, (lij) MX ■ ■ , a, (jk), (k(H)), (i5),...) + (m k). (3.14) 

The large angle soft terms can also be included in this structure. Two almost colour- 
connected unresolved partons are removed from the underlying colour ordering and we 
impose the same radiation pattern of these partons from the underlying colour ordering 
as in (3.14). The mapping applied to the primary antenna does not matter because this 
term only contributes in the soft limit where all mappings become identical. For six or 
more parton processes, the large angle soft terms can be generated by allowing the first 
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radiation to be between a pair of final-state partons [58]. In the special case of hadron- 
hadron initiated six parton scattering there will only ever be two hard final-state partons 
so the choice of final-state partons, between which the first unresolved parton is radiated, 
is unambiguous. However for five-parton processes there simply are not enough final-state 
partons and the first mapping should be of initial- final type. Four parton processes would 
require an initial-initial mapping, however in such a simple case there is no almost colour 
connected contribution. 

In terms of the radiation pattern described in (3.14) and Fig. 4, the first radiation 
takes place either between two final-state particles, or between one final-state and one 
initial-state particle, in the underlying ordering and the second radiation between the hard 
radiators in the antenna, j and /. If the hard partons being considered are both in the 
final-state then both unresolved partons are radiated between those hard partons. 

With the first mapping fixed to be of a particular type, the secondary antenna (which 
describes the radiation of the second unresolved parton) has common arguments across all 
terms, unlike in (3.14) where the arguments depend on the specific mapping inherited from 
the primary antenna. As a consequence, the secondary antenna can be factored out with 
the reduced matrix element so that, for example, in the final-final mapping case, the large 
angle soft can have the form, 

Y ■ Xl((ij), k, (iT)) M°(. . . , a, mk), (k(U)),b, ...), (3.15) 

where Y is a sum of large angle soft antennae. A similar structure emerges in the initial- 
final mapping case. The structure of this term can be understood with reference to Fig. 
4. For each region into which an unresolved parton can be radiated there is a term in Y 
given by, 



1 

2 



Sai/3 — SaIB 



(3.16) 



where / denotes the momentum of the unresolved parton i after the second mapping and 
for the example configuration considered in (3.14), 

(a,/3) = {(a,(g)),((5),(i0), ((7T),6)}, 

(A,B) = {(a,(^)),(^),(m)),({Hil)),b)}, (3.17) 

The pair of soft functions for each region comes with an overall sign depending on which 
region the primary antenna belongs to, i.e., a relative minus sign for the regions II and III. 
Collecting the terms in (3.14) and the large angle soft terms produces a block of terms for 
a given X®; the sum of such blocks constitutes da^ c NLO , 

d& NNLO = NnnLO J2 d $n+2(p3,---,Pn+4;Pl,P2) ^— Vj 

+ l - Xi(j, i, I) X$((ji), k, (zl)) MX . . , a, (MjX))A ■ ■ ■) 

- l - Xi(a, i,j) Xlm, k, I) MX . . , K), (Qjijk), (M),b, . . .) 
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\ X° 3 (l, i, b) Xl(j, k, {It)) M n °(. . . , a, (jk), (k(li)), (ib), . . .) 



( 5 (ij),*,(iO S ((ij)k)i(k(il))) 
~ ( S ai(ij) ~ S aiijij)k) ~ ( S bi(ii) ~ S bi({jX)k') 

Xl((ij), k, (71)) M° a (. . . , a, ((ijjk), (Mil)), b,...)\ 4 n \{ P } n ), (3.18) 



where the sum over i, k denotes the set of almost colour-connected pairs contained within 
the X® antennae in dcr^ L0 . The blocks in da^ c NLO have a common secondary antenna 
and so when integrated over the single unresolved phase space of the primary antenna, 
will stay together as a group of integrated antenna functions factoring onto a common 
unintegrated antenna. 

3.1.5 Colour disconnected subtraction term, da^^ LO 

For processes involving six or more coloured particles, colour disconnected double unre- 
solved configurations can arise where unresolved partons are separated by more than one 
hard radiator in the colour ordering. The single unresolved subtraction term da^^LO 
mits colour disconnected unresolved limits when the unresolved parton in the antenna is 
colour disconnected from any other parton in the reduced matrix element that can become 
unresolved. 

The subtraction terms in da^^ LO take into account all possible single unresolved 
partons and so for each colour disconnected pair of unresolved partons, j and m, where 
parton j lies in the antenna function and m in the reduced matrix element, there is the 
corresponding subtraction term where m lies in the antenna and j in the reduced matrix 
element. In this double unresolved limit both of these subtraction terms tend to the same 
value and da^ a NLO exactly double-counts the divergence of the matrix elements. 

To correct for this over-subtraction, a block of terms is introduced with the form, 

d& NNLO = -NnNLO Y\ d* n+2 (p3,...,p n +4;Pl,P2) 

Xi(i,j,k) X$(l,m,n) M° +2 (... ,I,K, -- ,L,N,---) J^(M„). 

(3.19) 



n+2 

X 



where the sum runs over pairs of partons separated by more than one hard parton in the 
ordering. 

S d 

The disconnected nature of the terms in d&^ NLO allows this block of terms to be 
analytically integrated over the two disconnected regions of single unresolved antenna phase 
space. The fact that both unresolved partons can be integrated out directly from the double 
real level dictates that this contribution should be reintroduced as part of the double virtual 
subtraction term. 
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3.2 Construction of the real-virtual subtraction term 

The real-virtual subtraction term must successfully remove the implicit IR divergent be- 
haviour and explicit poles of the one loop (n + 3)-parton matrix elements. In addition 
to imitating the physical one-loop matrix elements, the real-virtual subtraction term in- 
herits terms from the double real subtraction term after analytic integration over a single 
unresolved parton. da^ NLO contains three types of contribution: 

1. da T > a 

Terms of the type X® that cancel the explicit poles in the virtual one-loop (m + 3)- 
parton matrix elements. 

2. da T > b 

Terms that describe the single unresolved limits of the virtual one-loop (m+3)-parton 
matrix elements. 

3. da T > c 

Almost colour-connected contributions of the type ^3X3. 

The subtraction term for the real-virtual contribution can be further sub-divided into five 
contributions, 

d^NNLO = ^NNLO + d °NNLO + ^NNLO + ^NNLO + ^NNLO' (3.20) 

The structure and behaviour of each of these terms will be explained in the rest of this 
section, including whether the terms originate and terminate in the double real, real-virtual 
and double virtual contribution. 

3.2.1 Real-virtual mass factorisation term, ^fjNNLO 

The NNLO real-virtual mass factorisation contribution of Eq. (3.2) can be further divided, 

aa ij,NNLO — aa ij,NNLO + aa ij,NNLO' { 6 - zl -) 



such that da^^NLO is proportional to the (n + 3)-parton matrix elements, 



dx\ dx2 

X\ X2 



^f) C(e)T^ M (x 1 ,x 2 ) da^ NLO (x^ 1 H 1 ,x 2 ^H 2 ), 

(3.22) 

which will contribute to da^ NLO and da- N '^ LO , which contributes to da^ NLO . 

i*MF,l,b _ ,»AfF,l,t2 /ono\ 

aa ij,NNLO ~ aa ij,NNLO> {6.Z6) 

where, 

dxo 



x 2 
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(3.24) 

However, to simplify the book-keeping in the construction of the real-virtual subtraction 

MF,l,b 
ij,NNLO 



terms it proves useful to trivially rewrite daff^b b LO in the following way, 



da MF,l,b _ d& MF,lM A ~MF,l,b 2 Aa MF,l,b 3 ( ~ or.) 

aa ij,NNLO ~ aa ij,NNLO ^ aa ij,NNLO + aa ij,NNLO> {^.^0) 

where we have added and subtracted, 

AA MF,l,ln ((: rr t rr, _ f dx l dx 2 
d(7 ij,NNLO^lHl,&H 2 ) - j — — 

^ ^( e ) T ab]ab( X l^ X ^ d °ij,NL0( X lZl H l, X 2&H 2 ), 

(3.26) 



2tt 

M^nlo&Hi, &H 2 ) = -da^lo^H,,^). (3.27) 

The parton distributions that multiply these contributions to the partonic cross section 
are always fixed by the identity of partons i and j, which, in the initial- final (and initial- 
initial) configurations, are always of the same type as one (or both) partons involved in the 
antenna. However, in certain collinear limits, the type of parton can change, e.g. a gluon 
may transform into a quark or vice versa. Eqs. (3.26) and (3.27) describe these situations. 
In either case, the antenna is designed to correctly describe the singularity, but the parton 
emerging from the antenna (and therefore participating in the reduced matrix element) 
will be of a different type. In Eq. (3.26) a and b represent the species of the initial-state 
partons in the reduced matrix element in dcr^- NLO rather than the initial-state partons i and 
j that are involved in the antenna. In unresolved limits that do not change the identities 
of initial-state partons then a = i and b = j. However, in an identity changing initial-state 
collinear limit there will be a mismatch between the parton species in the antenna, i,j 
and those in the reduced matrix element, a, b. For example, consider a subtraction term 
d ®qg NLO th & t contains a term, 

d tf 9 ,NLO ~ ■ ■ ■ + A§(i„ 2 g , 3 ? ) M°(- %; ■■■) 4 n) ({p} n ), (3.28) 

where in the 2 9 ||3g limit, the initial state gluon transforms into an antiquark. In this case 
the initial-state partons involved in the antenna A®(l q ,2 g ,3q) are = {q,g) while the 
initial-state parton species in the reduced matrix element M^(- • • ; l q , 2q; ■ ■ ■ ) are (a, b) = 
(q,q). In other words, i = a = q, while the labels associated with parton 2 differ, j = g 
and b = q. In this example the mass factorisation kernel used in Eq. (3.26) is given by 

^ qq;qq{ X ^i x 2) m 

3.2.2 One-loop explicit pole subtraction term, d<x^'^ i0 

It was emphasised in section 3.1 that the construction of d<7^ LO follows the same lines as 
constructing ^cr% LO with a modified jet function and an additional parton. This matches 
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the interpretation of the mass factorisation contribution d<7 nnlo an< ^ so integrating the 
antennae in da^^ LO and combining with the mass factorisation kernels in d(7 nnlo wm 
generate precisely the same type of integrated antenna strings seen at NLO but with one 
additional parton, 

j-T,a __ f A ~S,a _ A ~MF,l,a 
aa NNLO ~ aa NNLO 17 N NLO 

ktRV f dxi dx2 ., 

= -N NNLO 2^ / d<f> n+ i(p 3 ,...,p n+3 ;x 1 p 1 ,x 2 p2) 

J X\ X2 



n+1 



~>n+l 



x jW 3 (l, . . . , n + 3) M° +3 (l, ...,„ + 3) 4 n+1) (Mn+i), (3.29) 

where M^ LO = C(e)M^§ L0 . It was shown in Sec. 2 that the poles of J^+2 are simply 
related to the poles of the (n + 2)-parton one-loop matrix elements. Using this fact, it is 
clear that the term d(T^-'^ i0 correctly subtracts the explicit poles of the real- virtual matrix 
elements. The only divergence remaining in the one- loop matrix elements is the implicit 
divergence arising from single unresolved configurations. 

3.2.3 Tree x loop subtraction term, d&^^ LO 

One-loop matrix elements factorise in implicit IR singular limits into two terms [68] which 
can be schematically understood as, 

1-loop — > (tree x loop) + (loop x tree). 

The first term is the product of a tree-level singular function factoring onto a one-loop 
reduced matrix element and is subtracted using a tree level antenna function and a one- 
loop reduced matrix element. This term by itself removes part of the implicit IR divergence 
of the one-loop matrix elements but introduces explicit poles associated with the one-loop 
(n + 2)-parton reduced matrix element. The explicit poles of this matrix element can be 
removed by introducing the appropriate integrated antenna string so that, after including 
the mass factorisation contribution da^^Q 1 , defined in Eq. (3.26), 



d(J NNLO = N NNLO 2^ / d$ n+ l(p 3 , • • • ,Pn+3] X 1 p 1 ,X 2 P2) ^~ 

J X\ X 2 On- 

n+1 

Xl(i,j,k)^5(l-x l )5(l-x 2 ) M* +2 (l,...,n + 2) 



x 

3 



+ cjJ%(1, . . . , n + 2) M° +2 (l, . . . , n + 2)| J^\{p} n ), (3.30) 

is free from explicit IR poles and where the sum is over the final-state unresolved partons. 
cj is a constant equal to unity unless n = and all particles are gluons in which case 
cj = 2 to account for both colour orderings of the gluons. This formula can be applied to 
final-final, initial-final and initial-initial configurations. 
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Al 


Al 




D\ 


D\ 


4 


4 


4 1 


pi 

^3 




G 3 


G 3 


G\ 


N x 


2 


1 


2 


3 


3 


2 


1 





3 


3 


2 


1 





M x 


1 


1 





2 


2 


2 





2 


2 


2 


2 





2 



Table 8: Number of colour-connected pairs Nx for the one-loop antenna X\, and the associated 
constant Mx- 



3.2.4 Loop x tree subtraction term, d<r^'^ i0 

The factorisation of one-loop matrix elements in IR limits also requires a term given by a 
one-loop universal singular function factoring onto a tree-level matrix element. To properly 
account for this contribution, a subtraction term is constructed from a one-loop antenna 
function and a tree- level reduced matrix element. The one-loop antenna function contains 
explicit IR poles which must be removed to ensure a finite total contribution. This goal is 
achieved by a difference of integrated antenna strings which contains the IR explicit poles 
of the one- loop antenna, and reflects the construction of [9] , 

y i r ..x o MUi,j,k) Q , . M\ (I, K) 

X s (x,j,k) =Si 3k ,iK 1 ^^-X^ 3 ,k) JiW ^ ) . (3.31) 

A block of terms can then be constructed to remove the remaining implicit IR diver- 
gence from the one- loop matrix elements, which incorporates the mass factorisation terms 
d&NNLO* an d ^nnlo 3 > defined in Eqs. (3.24) and (3.27) respectively, and is free from 
explicit poles, 

dx\dx2 .., si 



,.Tfe ktrv sr^ f dxi dx 2 . 

d(J NNLO = N NNLO 2^ / d$ n+ i(p 3 , . . . ,Pn+3]Xip 1 ,X2P2 

J X\ X2 



X 

3 



E 



n+l - — Dn+1 

X*(i,j, k)6(l - Xl )S(l - x 2 ) + J$>(i,j, k)Xl(i,j, k) 



M x X$(i,j,k)j£\l,K) 



M° +2 (.-.,I,K,...)4 n \{p} n ), (3.32) 



where Mx is a constant which depends on the species of one- loop antenna, X\. The string 
Jy) is constructed as a sum of Nx two-particle integrated antenna strings, 

N x 

J { x ] = E J 2 1] (hj), (3-33) 
(U)=i 

where Nx counts the number of colour-connected pairs of partons in the one-loop antenna. 
The values of Mx and Nx are displayed in Tab. 8. This subtraction term applies to 
all kinematic configurations and particle combinations. The mass factorisation kernels 
used to form the J x ^ come from da^NLO 2 an< ^ those used to generate come from 
d*^ nn'lo* ■ ^ ne term proportional to is derived from the integral of the double real 
subtraction term. Following the discussion in Sec. 2, the poles of a four-parton matrix 
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element integrated over the single unresolved phase space may be encapsulated by an 
integrated antenna string. Schematically, 



Voles 




= Voles 




= —Voles 


Ml 















(3.34) 

and therefore the single unresolved poles of a four-parton antenna are also described by, 
Voles J^Xl = Voles J?X^ = -Voles X\ . (3.35) 

In Sec. 3.1, da^f NLO was constructed to have no single unresolved limits such that, 



Voles 



[/ 



da 



~S,b 



Equivalently, 



Voles 



aa NNLO 



For each X^ or X% contributing to d<7^ LO , the associated block of iterated antennae in 
clcTyvATLO integrated over the single unresolved phase space therefore systematically gener- 
ates a singular contribution of the form, J^\i,j, k)X®(i,j, k), as required by Eq. (3.32). 

The one-loop antenna term, Xg, and the two-parton integrated antenna string, in 
Eq. (3.32) do not come from the double real subtraction term and so must be compensated 
for in the double virtual subtraction term. 

3.2.5 One-loop renormalisation subtraction term, d<x^'^ i0 

When introducing one-loop quantities it is important to ensure they are properly renor- 
malised to guarantee a complete cancellation of explicit poles. In the real-virtual subtrac- 
tion term two one-loop quantities are introduced: the one-loop reduced matrix elements 
in (3.30) and the one-loop antenna functions in (3.32). The one-loop matrix elements are 
renormalised at the renormalisation scale /U 2 whereas the one-loop antenna is renormalised 
at the mass scale of the antenna Sijk- To ensure both quantities are renormalised at the 
same scale, the replacement is made, 



NNLO 



-Voles 



da 



S,b 2 
NNLO 



(3.36) 



(3.37) 



X£(i,j,k) -> Xl(i,j,k) + ^ X$(i,j,k) 



Sijk | 



1 



(3.38) 



The systematic inclusion of these terms requires that each one loop antenna (whose pole 
structure is proportional to a tree-level antenna) used in the real-virtual subtraction term 
is compensated by a term proportional to (3q, such that a block of subtraction terms is 
generated of the form, 



r ]^T,b 3 ___ kI -rv \^ 
NNLO — JV NNLO 

n+1 



dx\ dx2 

Xi X 2 



d$n+l(P3, • • • ,Pn+3;Xip 1 ,X 2 P2) 
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x l3o\og(^—)xi(i,j,k)S(l-x 1 )S(l-x 2 ) 
j \\ s ijk\/ 

x M° +2 (... ,I,K,-..) 4 n \{p} n ), (3.39) 

where the sum is the same as in (3.32). The colour decomposition of /3q into 60 and 6o,JV F 
dictates to which orders in the colour decomposition these terms contribute, 

A> = b N + b , F N F , (3.40) 

where bo = 11/6 and 6o,f = —1/3. The terms in da^f^ LO originate in the real-virtual 
subtraction term and so, by integrating the three-parton antenna over the single-unresolved 
phase space, are reintroduced as part of the double virtual subtraction term. 

3.2.6 Integrated almost colour connected subtraction term, do-^J^LO 

The final term that contributes to the real- virtual subtraction term is mostly derived from 
the analytic integration of d&NNLO over the single unresolved antenna phase space with 
additional predictable terms to ensure an IR finite contribution. Integrating Eq. (3.18) 
over the single unresolved phase space and introducing three additional terms to ensure all 
explicit IR poles are cancelled yields, 

i-Tc K rRV sr^ f dxidx2 ,, , .If 
da NNLO = -Nnnlo >^ / d$ n+1 (p 3 , . . . ,p n+3 ;xipi,x 2 p2) e i 

oE f ( X °(. s ik) - X 3( s (ij)(jk))) 



2 

3 



{S(s ik , Sik, 1) — S(s^j^j k ^, Sjfc,^(ij)(jfe),jjfc)) 

~~ (S{ s aii s iki %ai,ik) ~ S{ s a(ij) , Sifc , 3 ; a(ij),ifc)) 

- (S(s k b,Sik,x kbyik ) -5(s (jfc ) 6 ,Sifc,a;( ifc ) Mfc ))^5(l - xi)S(l - x 2 ) j 

x X°(i,j,k) M° +2 (--- ,I,K,-..) 4")({p}n). (3-41) 

In this equation, the terms introduced to ensure IR finiteness are those involving the 
integrated antennae with mapped momenta, X®(suj\rj k \), X®(s a uj\) and X®(st k j\ b ). These 
terms must therefore appear in the double virtual subtraction term. For future reference 
we label these as, 

..Tci KrRV sr^ f dxidx 2 ,, , ,1 
da NNLO = -Nnnlo / d$ n+ i(p 3 , ...,p n+3 ;x 1 p 1 ,x 2 p 2 ) ~ 

~^ J Xl X 2 O n+ l 
9 ^2 [^3°( s (ij)(jfc)) + X i( s a{ij)) + X i( s {kj)b)] 



X 

2 

3 



x X$(i,j,k) M n ° +2 (- ,I,K,---) 4 n) ({p}n), (3-42) 
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,.Tm , rRV ST^ f dxi dx2 ,, , \ 1 

da NNLO = +N N NLO / d$ n+ i(p3, ■ ■ ■ ,Pn+3;XiPi,X2P2) ~ 

~{ J Xi X 2 O n+ i 

x E -*3(*(«)0-*)) *3°(^\*0 M„° +2 (... ,/,#,-■■) jW({p}n). (3-43) 



The integrated soft function for the final-final mapping is denoted by S(s ac , Sj&, x a c,jA;)- An 
explicit expression can be found in Ref. [58]. As discussed in Sec. 3.1.4, for five-parton 
processes it is necessary to use a soft function with an initial-final map. In this case one 
finds an analogous equation to Eq. (3.41) involving the integrated initial-final soft function, 
S IF (s ac ,s I K,yac,iK), which is given explicitly in App. B. 

3.3 Double virtual subtraction term structure 

The double virtual contribution to the pp — > n-jet cross section involves the two- loop 
(n + 2)-parton matrix elements which have no implicit IR divergence in any regions of the 
appropriate n-parton phase space. All that remains is to reintroduce the integrated forms 
of the appropriate subtraction terms, and combine them with the double virtual mass 
factorisation terms, so that the explicit IR poles of the two-loop contribution are cancelled. 
The main result of this section is that the double virtual subtraction term is constructed 
from three terms, 

^ij,NNLO = ^ij'^NNLO + d ^fj',NNLO + d<5 ij, N NLO ■ (3AA) 

The following sections will demonstrate the origin of each of these terms. This construc- 
tion reflects the dipole-like singularity structure for two-loop amplitudes made apparent in 
Catani's two-loop factorisation formula [64], 

Voles{M 2 n (l,... ,n)) = 2/«(e; 1, ■ ■ ■ ,n)^(l,-- - , n )-^M°(l-- - , n) 
- 2lW(e;l,--- ,n) 2 M°(l,-- - ,n) 

+ 2e ~ £J W^i (t +*)4 1) (26;l - ,n)M°(l... ,n) 

+ 2H^(e)M^l,... ,n), (3.45) 

with I { n\e) given by Eq. (2.34) and where the form of the hard function, H^ 2 \e) and the 
constant K depends on the particle content and order in N under consideration. 

3.3.1 Double virtual mass factorisation terms, d^/yfjvxo 

The general form of the double virtual mass factorisation contribution was given in Eq. (3.3). 

(2) 

To simplify the construction of the double virtual subtraction term, F^' kl (zi, z 2 ) may be 
decomposed in the following way, 

rgk(*i, *>) = rgk(*i, z 2 ) - ^ kl ( Zl , Z2 ) + \ [rgj", ® r% t ] ( Zl , z 2 ), (3.46) 



such that, 



rgiz^i, z 2 ) = Tf k \ Zl )5 3l 5(l - z 2 ) + Tf(z 2 )5 lk 5{l - Zl ), (3.47) 
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where, 



1 

2e 



p^z) + ^p%{z) ) . 



The double virtual mass factorisation contribution can be recast into three terms, 

,»MF,2 _ A ~MF,2,A ,.MF,2,B A ~MF,C 
a<J ij,NNLO ~ aa ij,NNLO + aa ij,NNLO + aa ij,NNLO> 

where the individual contributions are given by: 

'a,N 



aa ij,NNLO 



I" dzidz2 ( 

J Zi Z 2 V 



do. 



MF,2,B 
ij,NNLO 



+ 



d<7 



MF,2,C 
ij,NNLO 



dz\ dZ2 
Z\ Z 2 

dz\ dz2 
Z\ z 2 
dz\ dz 2 
z\ z 2 
dz 2 

Z-2 



2vr 
a s N 
2ir 

2ir 



C(e)r%(da 
C(e)T^ kl da T 



T kl,NLO 



— —aaki,LO 



(3.48) 



(3.49) 



(3.50) 



ij;kl^ u kl,NLO 



(i) 



[d Zl dz 2 (a s N\ =., 2 =(2) , 



(3.51) 
(3.52) 



and the convolution of two functions f(x\,x 2 ) and g(yi,y 2 ) is defined in the following way, 
[/ ® g] (zi,z 2 ) = J dx 1 dx 2 dy 1 dy 2 f(x 1 , x 2 )g(y 1 ,y 2 )5(z 1 - x 1 y 1 )8(z 2 - x 2 y 2 ). (3.53) 
Each of these terms naturally fits into the corresponding piece of the double virtual sub- 



traction term do 



U,A 



' ij,NNLO' dvfj^?NLCn and ^ij%NLO^ to ren der each contribution free from 
initial-state collinear poles. 

3.3.2 d<xJ^ LO subtraction term 

The first double virtual subtraction term is generated by integrating two terms from the 
real- virtual subtraction term and combining these with the appropriate mass factorisation 
term. The first term to be integrated comes from d<r^'^ i0 , defined in Eq. (3.30) as the 
collection of terms proportional to the one-loop matrix element. After integration, this 
contribution yields a sum of integrated antenna functions factoring onto a one-loop matrix 
element. The second term comes from the integrated form of da^/jy LO , 



~U,C 



aa ij,NNLO — ~ JV NNLO 



{iJ} 



Bo 



M 2 



d<S> n {p3, . . . ,Pn +2 ; ZlPl, Z 2 p 2 ) — 

zi z 2 S n 



-1) Xl{ Sl] )M» n+2 {l r .. ,n + 2)jH({p}, 



(3.54) 



where M^Yilo = ^( e )-^NNLO = ^( e ) 2 -^NNLO- ^he bracket may be expanded to produce 



The 



two terms. The first term (proportional to (|sjj| e )) will become part of do 
second, proportional to —1, is combined with the one-loop matrix element term from the 



U,C 

NNLO- 
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integral of dcr^ LO and the mass factorisation term du^^LOi to produce the first double 
virtual subtraction term, 

da ij'jVNLO = -^NNLO Yl l ~~ d ^n(P3, • • • ,Pn+2] Z-tfH, Z 2 p 2 ) ^~ 

J z l z 2 <->n 

x j£ 2 (l, - ,n + 2)^ +2 (l,... ,n + 2)-^ M n ° +2 (1, • • ■ , n + 2)) 

(3.55) 

The poles of the integrated antenna string are directly related to the poles of Catani's 
one-loop insertion operator, therefore this contribution to the double virtual cross section 
contains the 1/e 4 and 1/e 3 singular contributions given in the first line of Eq. (3.45). 

It is important to note that da^^ LO does not contain precisely the same singularities 
at 1/e 2 and above as the first line of Eq. (3.45) because the finite difference between 
and IrP is formally of order e°. These differences ultimately cancel against similar terms 

in d& NNLO and d °NNLO- 

3.3.3 d&j^ NLO subtraction term 

The second double virtual subtraction term is formed from the integral of the remaining 
term from da^ ^ NLO (3.30) proportional to M", the integral oida^^ LO (3.43) and da^ NLO 
(3.19), and the mass factorisation contribution da NN ' L ' . The resultant subtraction term 
is given by, 



da ij'NNLO = -Nnnlo Yl [ —— d$ n (p 3 , . . . ,p n+2 ; z 1 p 1 ,z 2 p 2 ) 

n •* 1 2 

x "[4+ 2 (V-- ,n + 2)® J%(1,--- ,n + 2)]( Zl ,z 2 ) 



1 



2 

x M° +2 (l,--- ,n + 2) 4 n \{p} n ). (3.56) 

It is easily seen that this expression contains the pole structure of the second line of 
Eq. (3.45) with additional singular contributions arising from the finite differences between 
Jn^ and 1$ at O(e ) and higher orders. 

Note that the combination of Eqs. (3.55) and (3.56) reproduce the pole structure of 
the first two lines of Eq. (3.45) up to 0(l/e). 

3.3.4 da-pJ NLO subtraction term 



The third double virtual subtraction term is generated by integrating appropriate terms 

T NNLO (3-ll)i do^^LO 



in the remaining double real and real- virtual subtraction terms, da^ h ^ LO (3.11), d<rj' bl 



(3.32), da^^LO (3-43), the term proportional to |sjj| _€ generated by expanding the bracket 
in Sec. 3.3.2 and the mass factorisation contribution da^^LO (3-52). The final subtraction 
term is given by, 

da NNLO = -NnNLO ^2 f —— d$> n (p 3 , . . . ,p n+2 ; ZlPl, Z 2 p 2 ) ^- 
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Ji 2 ] 2 (l,-- - ,n + 2) M° +2 (l,--- ,n + 2) 4 n) ({p}n), (3.57) 

(2) 

where we have introduced the double unresolved integrated antenna string, J n + 2 as a sum 

(2) 

over double unresolved integrated dipoles, J \ , 

J%(l,---,n + 2) = y £j?\hj)i (3-58) 

(id) 

and the sum runs over colour connected pairs of partons in the (n + 2)-parton ordering. 
The double unresolved integrated dipoles depend on the type of particle in the scattering 
process, the order in the colour decomposition under consideration and the kinematic 
configuration (FF, IF or II) of the dipole. Schematically, the final-final double unresolved 
integrated dipole has the form, 

jf } (I, J) = c FF X°( Sl j) + c FF x!(s u ) + c FF Xl{ SlJ ) 

+ c FpP -f (i^l) ^XKsu) + c FF Xl(s LJ ) *s (ay), (3-59) 

where c FF , are constants associated with each integrated dipole. For an initial-final in- 
tegrated dipole, the form of the integrated dipole is modified by the use of initial-final 
integrated antennae and the inclusion of the necessary mass factorisation contribution to 
remove all initial-state collinear poles, 

J? (1,1) = c{ F xZ(s- u ) + cl F xZ(s u ) + 4 F Xl(s^ 

+ c IF Xl{s- u ) X°(s- U ) ~ T%\zi)6(l - z 2 ), (3.60) 

where i labels the species of parton in the initial-state and k is the parton involved in the 
matrix element in Eq. (3.57). The initial-initial integrated dipole is formed from initial- 
initial antennae and includes a mass factorisation contribution with non-trivial z\ and Z2 
dependence, 

jf (U) = ^m+^m+^^ + ^M)^^) 

+ <%xl<,m) ® ^3(^12) - r%( Zl )5(i - z 2 ), (3.6i) 

—12) 

where the i,j labels carried by - w are the species of initial-state parton carried by the 
cross section and k, I denote those carried by the matrix element in Eq. (3.57). 

3.3.5 The full double virtual subtraction term, ^YiNLO 

The full double virtual subtraction term is given by the sum of dcr^ LO (3.55), da^ LO 
(3.56) and da^ LO (3.57) so that, 

= --^NNLO ^2 f —— d$„ (p 3 , . . . , p n+2 J Ztfi , Z2P2 ) { 
„ J Z\ Z2 Jn { 
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Figure 5: NNLO double virtual structure. Integration of the tree-level double real emission over 
the double unresolved phase space, plus the integration of the one-loop real-virtual contribution 
over the single unresolved phase space, generates the double virtual singularity structure described 
in Eq. (3.62), for each colour connected dipole pair /, L. The one-loop diagram on the right of this 
figure implicitly includes the /3q term in Eq. (3.62). 

J%(1,--- ,n + 2)(iVC(l,- ,n + 2)-^ M n ° +2 (1,--- ,n + 2)j 
+ \ [JnUh ■ ■ ■ . n + 2) <g> j£> 2 (l, • • • , n + 2)] {z u z 2 ) M° +2 (l, • • • , n + 2) 
+Ji 2 ] 2 (l,--- ,n + 2) M n ° +2 (1,-.. ,n + 2)\4 n \{p}n). (3.62) 

This master equation encapsulates the IR poles generated by the integration of the tree- 
level double real emission over the double unresolved phase space, plus the integration of the 
one-loop real-virtual contribution over the single unresolved phase space and is summarised 
in Fig. 5. 

As previously discussed, the first term in Eq. (3.62) reflects the pole structure of the 
first line in Eq. (3.45). The second term similarly reflects the structure of the second line 
in Eq. (3.45). The third term introduces the double unresolved integrated antenna string, 

(2) 

which is summed over colour-connected dipoles. This term corresponds to the terms 
proportional to J^. 2 (2e) and H^- 2 \e) in the last two lines of Eq. (3.45). 

Of course, it is known that for six or more partons, there are potentially multiple par- 
ticle correlators in the 1/e contribution. Nevertheless, once the scale dependent functions 
proportional to |sij| _e are expanded, the difference between Eq. (3.62) and the two-loop 
matrix elements is finite. 

Note also that when adding together all contributions in the double-virtual chan- 
nel, we observe that contributions of order e or higher to the one-loop amplitude can- 
cel, even though they are multiplied with divergent factors in individual terms (one-loop 
self-interference and integrated real-virtual subtraction terms). This cancellation can be 
understood in detail from the NNLO infrared singularity structure [64,69]. 

4. Quark-antiquark scattering to gluonic dijets at NLO 

In this section we review the IR structure of gluonic dijet production via quark-antiquark 
scattering at NLO (to all orders in ./V and setting Np = 0) using the antenna subtraction 
method. 

Here we first review our notation for amplitudes involving quarks. In Sees. 2 and 3 the 
symbol M l n was used to denote a generic /-loop, n-parton, squared partial amplitude. In 
previous works [58,59,67] the symbol A l n was used to denote the /-loop, n-parton all-gluon 
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squared partial amplitudes. Following this example, squared partial amplitudes involving 
one quark- antiquark pair and (n — 2) gluons will be represented by the symbol B l n and we 
will reserve C l n (D l n ) for processes involving two non-identical (identical) quark-antiquark 
pairs and (n — 4) gluons respectively. 

In our notation, the leading colour /-loop amplitude for (to + 2)-parton scattering with 
one quark-antiquark pair in the initial-state is given by, 



°m+2 



(te,A !lM , J }) = 2»/y(^) i £ 

V 1 J aes n 



x ( T Mi)...rv™)).. B l m+2 (l q ,a(3) g ,--- ,a(m + 2) g ,2 q ), (4.1) 

where the permutation sum S m is the group of permutations of m symbols. The SU(N) 
algebra is normalised according to the convention, Tr (T a T*) = 5 ab /2. Squaring and sum- 
ming over helicities and colours gives the leading colour Z-loop qq — >■ m gluon contribution 
to the n-jet cross section, 

da = A/"4 +2 d$ m (p3, • • ■ ,Pm+2-,Pl,P2)^ 

x E ^ + 2(i 9 ^(3)„...,cT(m + 2) 9 ,2 9 -)j( m )(M m ), (4.2) 

where B l m+2 denotes the /-loop, (m + 2)-parton, colour-ordered helicity summed, squared 
partial amplitude. The tree-level, one-loop and two-loop squared partial amplitudes are 
defined according to: 

B° m+2 (a) = \Bl +2 {o)\\ 

B 1 m+2 (a) = 2Re[B°J +2 (a)B 1 m+2 (a)], 

B 2 m+2 {a) = 2Re[B° m \ 2 (a)B 2 m+2 (a)] + |^ +2 (^)| 2 , (4.3) 

where a denotes a given colour ordering. The normalisation factor M} n+2 includes the 
average over initial spins and colours and is given by, 



AC« = Af io X ( ^-j (4.4) 



where for the qq — > gg process, 

The coupling g 2 has been converted into a s using the factors C(e) and C(e), 

9 2 NC(e) = C(e). (4.6) 

For low multiplicity matrix elements (m ^ 3), the sub- leading colour contributions can be 
written as an incoherent sum of squared partial amplitudes. 
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4.1 Leading order cross section, d&Lo 

The leading order contribution to qq — > gg scattering is given by, 

&Vqq,LO = MlO J &$2 (P3 , P4\ Pi , P2) y 

' P(*J) 
0/ 



Bt(i q ,ig,j g ,2 q ) - ^ B^i q ,i g , jg ,2 q )\ 4 2) (P3,P4 



(4.7) 



where P(i,j) denotes the two permutations of the gluons 3 and 4 in the colour ordering. 
The sub-leading colour contribution is formed from the square of a coherent sum of colour- 
ordered matrix elements, 

B° 4 (i q ,i g ,j g ,2 q ) = \B° 4 (i q ,i g ,j g ,2 q ) +B° 4 (i q ,j g ,i g ,2 q )\ 2 (4.8) 

Here the three gluon coupling involving gluons i and j drops out and this contribution is 
sometimes referred to as QED-like, or as involving Abelian gluons. 

4.2 Virtual cross section, d<x]^ LO 

The Np independent one-loop contribution to the qq — > gg subprocess is given by, 

d &qq,NLO = Ml J d$2(P3,P4;Pl,P2)^^<*(l - Xi)S(l - X 2 ) \ ^ 



1 



B 4 {lq, ig,jg, 2 q ) ^ 2 



where the overall factor is, 

Ml = M LO 



1 



P(i,j) 



+ ^4 B l(l q ,i g ,j g ,2q) \ 4 2) (Pi>Pj) 



2ir 



C(e). 



(4.9) 



(4.10) 



The various contributions are formed from the projection of one- loop partial amplitudes 
onto tree-level amplitudes. The first sub-leading colour contribution is split into two terms, 
B\ a and B\ b to reflect that fact that the explicit poles of these contributions factor onto 
different tree-level matrix elements. The poles of these contributions may be expressed in 
terms of integrated antenna strings, 



Voles Bl(l q ,ig,j g ,2q) 
Voles B\ a {l q ,i g ,j g ,2 q ) 
Voles Bl b (l q ,ig,jg,2 q ) 



J 4 ^ (lg> igtjgi ^<?) B 4 (l q , ig,jg, 2 q ), 
J 2 (lg>^<z) B 4 (l(j, ig, jg, 2q)i 



(«^3 0-q> ^<?) + J 3 (iqijgi ^q)) B2(i-q,ig,j g ,2 q ) 
+ji\lqi%) Bl(l q ,ig,jg,2 q ), 
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Voles 



Bl(lq, ig,jg, 2 q ) 



^2 \^q^q) B*l{l q ,ig, jg,2 q ). 



(4.11) 



Following the discussion of Sec. 2, the integrated antenna strings are formed from sums of 
integrated dipoles, 

J 4 \^-qi igijgi %) = ^2 \^-qi h) + ^2 \^gijg) + *^2 \jg-> 2g)i 
^0-qi igi %) = J 2 ^0-qi h) + *^2 \jg-> 2g); 



\jg^q) ~ *^2 ^ Jfir): 



(4.12) 



where the relevant integrated dipoles are given in Tabs. 4 and 5. 

Using these relations, the poles of the Np independent virtual cross section from the 
channel qq —> gg can be rewritten in the form, 



Poles (da^ NLO ) = -Ml j d$2(P3,P4;Pi,P2) 



dxi dx2 1 

X\ %2 2 



J 4 \^-qi^gi jgi^q) 0-q^gj jg, 2<j) 



1 

( N 2 + l 



J 2 ' '(^9' 2 9 ) B^(lg,ig, jg, 2g) + J3 ^(Iq, Ig, 2g) Bl(l q ,ig, jg, 2g) 



J^O-q^q) ^4(lq)3 fl ,4 fl ,2g) 



(4.13) 



4.3 Real cross section, d<r^ 



NLO 



The real emisson contribution to the Np independent qq — > gg NLO cross section comes 
from the tree-level five-parton qq — > ggg process, 



d °qq,NLO = ^5 d ^3(P3,P4,P5-,Pl,P2) ^ Yj 



B5(lq,ig,jg,k g ,2 q ) 

P(i,3,k) ' 

^(^q^gdg^g, 2g) + — f j-j^ ^ B®(i q ,ig,jg,kg,2q) 



^2 3) (Pi,Pj,Pfe)|, 

(4.14) 



where P(i,j, k) is the set of six permutations of the gluons in the colour ordering and the 
overall factor is given by, 



\r0 K r fcx s N\C(e) 

=Mlo {-2^)wy 

The sub-leading colour matrix element B®(l q ,ig,jg,kg,2 q ) is given by, 



(4.15) 



Bs(iq,ig, jg, kg,2 q ) — 

\B 5 (l q , ig, jg, kg, 2 q ) + B^jg, ig, kg, 2 q ) + B^jg, kg,ig, 2 q )\ 2 , (4.16) 
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such that the gluon i behaves in an Abelian fashion. The most sub-leading contribution is 
given by the QED-like matrix element, formed by averaging over all colour-ordered matrix 
elements, 

2 



^5(1?' *S> Jfl' kg, 2g 



^(iq, i 9 , j 9 , 2^) 



(4.17) 



All gluons in this contribution are QED-like and are only colour connected to the quark- 
antiquark endpoints. 

4.4 Real radiation subtraction term, da^ LO 

The subtraction term for the real cross section presented in Eq. (4.14) is given by, 



d °qq,NLO = ^5 d $3 (P3 , PA, P5 ] Pi , P2) ^ ^ 



1 

iV2 



+ 7 



1 /A 2 + 1 



3! \ A?" 4 



3! 

P(i,j,k) 

4(l,i,j) B2(l q ,(if) g ,k g ,2g) jf\ P(ij) ,P k ) 

+ fl(i,j,k) B° 4 (l q ,{ij) g ,(jk) g ,2 q ) 4 2) (p m ,p m ) 
+ 4(2,k,j) B° 4 (i q ,i g ,(jk) g ,%) 4 2) (pi,p (jk) ) 

A° 3 (l,i,2) B° 4 (t q ,j g ,k g ,%) 4 2 \ Pj ,p k ) 
+ 4(1, j,k) B%(l q ,i g ,(jk) g ,2 q ) 4 2) (pi,p W ) 
+ 4(2,k,j) B° 4 (i q ,i g ,(jk) g ,%) 4 2 \ Pi , P(jk) ) 

A° 3 (l,i,2) B° 4 (l q ,j g ,k g ,2 q ) J {2) (p 3 ,p k ) 
+ A° 3 (l,j,2) B° 4 (t q ,ig,~k g ,2 q ) J [2 \pi,p k ) 
+ A° 3 (l,k,2) B° 4 (i q ,i g ,j g ,%) 4 2) (PhPj 



(4.18) 

4.5 Virtual subtraction term, d(T^ LO 

Integrating the real subtraction term over the single unresolved phase space and combining 
with the appropriate NLO mass factorisation kernels allows the virtual subtraction term 
to be constructed from integrated antenna strings, 

dx 2 



d°qq,NLO = 



I d$ 2 (p3,P4;Pl,P2) J ~xt 



P(iJ) 



X2 



J 4 \lqi ig-> 3gi 2<j) B^O-qi ig-iJg-i 2<?) 



1 

'A 2 " 



\^qi 2g) ^4(1<J) 2g) + ^3 ''(Igj 2 q ) B 4 (l q , i g ,j g , 2 q 



'A 2 + l 
2A 4 



*^2 ^(^9) ^g) B 4 (l q , lg,jg, 2 q ) 



4 2) (Pi>Pj)> 



(4.19) 
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From the definition of these strings the connection to the unintegrated subtraction term in 
Eq. (4.18) is clear. This expression exactly matches that of Eq. (4.13), i.e., 

Voles (da^ NLO ) - Voles (da^ NLO ) = 0. (4.20) 

5. Quark-antiquark scattering to dijets at NNLO at leading colour 

The leading colour NNLO contribution to dijet production has previously been studied 
using the antenna subtraction formalism for purely gluonic channels at the double real [67], 
real-virtual [58] and double virtual [59] levels. 

This section will focus on the leading colour NNLO corrections to dijet production from 
quark-antiquark scattering. These corrections include the double real tree-level contribu- 
tion qq gggg, the real-virtual one-loop contribution qq — > ggg and the double virtual 
two loop contribution qq — > gg, all of which will be examined using the methods of Sec. 3. 

5.1 Double real cross section, d<x^y i0 

The leading colour contribution to the cross section from the qq — > gggg subprocess can be 
written in the form, 

dtqq R NNLO = K d$ 4 (p 3 , ■ ■ • , P6', Pi , P2) J { 

' P(i,j,k,l) 

X B%(i q ,i g ,jg,kg,lg,2q) (Pi , Pj , Pk , Pi) , (5-1) 

where P(i,j, k, I) denotes the 24 possible permutations of gluons in the colour ordering and 
the overall factor is given by, 

^- Mlo \-w)cw (5 - 2) 

It is convenient to rearrange the 24 squared amplitudes present in Eq. (5.1) into three 
terms, 

M(i(iq,ig,jg,kg,lg,2q-) = Xg (l q , 3 g , 4 g , 5 g , 6 g , 2q) 

P(i,j,k,l) 



+ ^6 (!<?> 3g, 5 g , 4g, 6g, 2 



+ X$(i q ,3 g ,4 g ,6 g ,5 g ,2g), (5.3) 

where each X® contains eight colour ordered squared amplitudes obtained by the four cyclic 
permutations of the final state gluons and their line reversals: 

"^6 (!<?' 4 fl , 5g, 6g, 2q) = B®(i q , 3g, 4g, 5g, 6g, 2g) + Bq (l q , 6g, 5g,4g, 3g, 2g) 

+ -^6 (lg> 4 g , 6 9 , 3 g , 2g) + Bq (lq, 3 g , 6 g , 5 g , 4 g , 2q) 

+ B$(l q , 5g, 6g, 3g, 4g, 2q) + 4g, 3g, 6g, 5g, 2q) 

+ B° 6 (i q , 6 g , 3 g , 4 g , 5 g , 2 q) + B° ( 1 „ 5 g , 4 g , 3 g , 6 g , 2 q ). (5.4) 

It is sufficient to apply the subtraction technique to one block of orderings, the other two 
are related by symmetry and contribute numerically the same result after integration over 
the phase space. 
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5.2 Real-virtual cross section, d(T^^ i0 

The leading colour one-loop contribution to the real-virtual cross section for the qq —> ggg 
subprocess is given by, 

&°qq,NNLO =•^5 / d $2 (P3 , PA, Ph\ Pi , P2) }^ 

J 6. X\ X 2 

x B b (\ q ,i g ,j g ,kg,2q)8(l - x{)8(\ - x 2 ) J 2 (Pi,Pj,Pk)- (5.5) 

where P(i,j,k) denotes the set of six permutations of gluons in the colour ordering and 
the overall factor can be derived from (4.4), 

,.j Kr /a s N\ 2 C(e) 2 . . 

and 

MZlo = Mlo (^) 2 ^ = MZ LO C(e). (5.7) 

5.3 Double virtual cross section, ^Y^nlo 

The leading colour contribution to the two-loop qq — > gg cross section is given by, 



d °qq V NNLO =M% I d* 2 (P3,P4;Pl,P2) h^T^T Yl 

J 1 Z2 P(i,j) 

x Bl(i q ,ig,j g ,2g)S(l - xi)S(l - x 2 ) 4 2) (Pi,Pj), (5.8) 

where the overall factor is, 

Ml=N LO {*£-) 2 C{e)\ (5.9) 

As in [58,59] and Eq. (4.3), the two-loop matrix element contains the projection of the 
two-loop amplitude onto the tree-level amplitude and the self-interference of the one-loop 
amplitude. 

5.4 Double real subtraction term, d&^ NLO 

The leading colour contribution to the squared matrix element is an incoherent sum of 
squared colour-ordered partial amplitudes. The sum over colour orderings can be parti- 
tioned into three blocks of orderings as described in Sec. 5.1, such that a subtraction term 
may be constructed for a block of orderings, rather than the entire squared matrix ele- 
ment. Following the general discussion of Sec. 3, the NNLO subtraction term for the block 
of orderings in Eq. (5.4) is given by, 



^NN 6 LO =Ar 6^(P3,..-,P6;Pl,P2) ^ £ j 

' P c (ijkl) I 
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+fi(hj,k)B^(l q , (h3)g,ti k )g, l g,2q)J2 3) (P(ij)iP(jk),Pl) 

+/ 3 °(j, k, Z) B%(l q ,i g , 2g) 4 3 \p(i),P( jk ),P(kl)) 

i, j) /%, Z fl , 2g) 4 3 \p(ij),Pk,Pi) 

+4(2, I, k) B°(l q ,i g ,j g , (kf)g, %) 4 3) (^Pj'P(kl)) 

+/s (*> k, j) B%(l q , (lk) g , (kj) g , i g , 2 q ) jf } (p {lk) , , K ) 
+fl{k,j,i)B^(l q ,lg, (kj) g , (ji) g ,2g) 4 3 \pi,p {kj) ,p {ji) ) 
+4(1,1, k)B%(l q , (lk) g ,jg,ig,2 q ) jf\p(ik),Pj,Pi) 
+d§(2, i, j) Z g , 2g) (pi,p k ,P{ji)) 

+ F 4, a (h3i k , 5 4(ig> (ij^g, {jkl) g , 2 q ) jf ] {P(ijk),P(jkl)) 

+ F 4,b(h3> h k ) B l(l q , (ijl)g, U l k) g , 2q) J2 ] {P(ijk),P(jlk)) 

+Bl a (l, k,j, i) B2(i„ 2 5 ) 4 2) {P(ikj),P(kji)) 

+F 4 ° 6 (Z, fc,i, j) B2(i„ (ifci) p , (kij) g , 2 q ) J 2 (2 ] (p {m) ,p {kij) ) 
+Dl(l,i,j, k) B%(l q , (ijk) g ,l g , 2 q ) 4 2) ( P{tjk) ,pi) 
+Dl(2, I, k,j) B%(i q , i g , <jkj) g , %) 4 2) (PhP(lkj)) 
-A%1, h k, 2) B° 4 (l q ,j g , Ig, %) 4 2) (Pj,Pi) 

-fKi,j, k) fim, (jk),i) Bt(i q , (mjk)) g , (mi) g , 2 q ) jf\p mm) ,p {{lk)l) ) 

-f!(j, k, i) fi(i, (jk), (id)) b ^ (7(jk)) g , ((Jmi)) g , 2 q ) 4 2) (p mW ),P( W i)) 

- fi(i,j, k) fim, I, (jk)) Bl(l q , Uifftg, ®) s , 4 2) (p m t),p iWl) ) 

-fl(i,k,j)fim,(^),i)Bt(i q ,((im 

-f$(k,j, i) fi(l, (kj), (ji)) Bl(l q , (7(kj)) g , ((kj)(ji)) g , 2 q ) 4 2) (P(k J ),P((k J m)) 

- f$(l, k,j) f!((lk),i, (kj)) Bl(l q , ((lk)i) g , (&jji) g , 2 q ) J?\p mi) ,p [{kj)i) ) 

-4(1, i, j) D° 3 (l, (if), k) B%(i q , ((ijjk)g, l g , 2 q ) 4 2) (P((ij),k),Pl) 
-f$(i,j, k) D$(l, (ij), (jk)) B° 4 (I q , ((ijj(Jk)) g , l g , 2 q ) 4 2) (Pm(jk)),Pl) 
-4(1, k,j) D° 3 (l, (kj), i) Bl(l q , ((kj)i) g , l g , 2 q ) 4 2) (Pm)i),Pl) 

-4(2, 1, k) D° 3 (2, (lk),j) Bl(l q , i g , m)j) g , %) 4 2) (P((lk) j ),Pi) 

-fl(l, k,j) D° 3 (2, (Ik), (kj)) Bl(l q , ig, ((lkj(kj)) g , %) J?( Pl ,p m{m ) 
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-4(2, j, k) Dl{% (jk),l) B° 4 (l q , i g , ((jk)l) g , %) 4 2) (pi,p mt) ) 

+ Ajfti, i, 2) A° 3 (t, k, % B° 4 (t q ,] g , l g , %) jf\p~ pPl ) 
+ A° 3 (l, k, 2) A° 3 (l, i, % B° 4 (l q ,] g , l g , %) 4 2) {p,, P! ) 

+\ fi(i, i, k) fim,j, (jk)) B° 4 (i q , {(u)j) g , mjj) g , 2 q ) 4 2 \p((u)j)>pm)i)) 

-\ 4(1, I, i) f 3 ((U),j, k) B° 4 (t q , (5oJj fl) (jk) g , 2 q ) 4 2 \p mj) ,p Uk) ) 

-\ 4(2, i, k) fl(i,j, (ki)) B° 4 (i q , (ij) g , (j(kt)) g , %) 4 2) (p<ij),pum) 

{S(u)i(ik) - S((ii)j)i(j(ik))) - {Sii(ii) - s ii((u)j)) - (S2i(ik) - S 2 i(j(ik))) 

XfV((U),j,(lk)) Bl(l q ,(^)g,(j(lk))g,2 q ) ^(p^PWk))) 

+\ fi(i, ij) Jim, k, m Bi(i q , mk) g , mk) g , %) j {2 \p mk) , vmk) ) 



-- 4(1 i, fim, k,j) Bl(l q , ((il)k) g , (kj) g , 2 q ) 4 \p((il)^P{kj)) 

-\ 4(2, ij) fi(i, k, m B° 4 (i q , (Jk) g , (Hij)) g , %) 4 2) (P(ik),P(km) 

(S(li)i(ij) - S((U)k)i{k{ij))) - (Sii(u) - £li((Ji)fc)) - (S 2 i(ij) - S 2 i(k(ij))) 

xfl((li), k, (ij)) Bl(l q , ((li)k) g , (k(ij)) g , 2 q ) 4 2) (P((ii)k),P(k(ij))) 

+ l - 4(1, k,j) 4(1, i, (kj)) Bl(l q , fikjji) g , lg, 2 q ) jf ] (p m) i),Pl) 

-\ fiU, k, I) 4(1 i, W) B° 4 (l q , (7(jkj) g , (H) g , 2 q ) 4 2) (P(iUk)),P(ki)) 
-\ A%(1, k,2) 4(l,i,j) B° 4 (l q , (jj) g j g , %) 4 2) (P(i,j),Pi) 

(Slk(jk) - Slk(i(jk))) ~ (S(jk)k(kl) - S(i(jk))k(kl)) - (Slk2 - Si k2 ) 

x4(l,i,(jk)) Bi(l q ,(7(j^) g ,(kl) g ,2 q ) J { 2 2 \p {l(jk)) ,p [kl) ) 

+ l - 4(1, %, j) 4(1, k, m B° 4 (t q , mk) g , i g , 2 q ) j { 2 2 \p mk) ,pi) 

-\ fi(j, i, I) 4(1 k, (ji)) B° 4 (l q , (k(ji)) g , (U) g , 2 q ) 4 2) (P(k, {j ,i)),P(i,i)) 
-\ A%(l,i, 2) 4(1, k, j) Bl(l q , (kj) g , l g , %) 4 2) (P(kj),Pf) 



(Su(ji) - Su(k(ji))) ~ ( s (k(ji))i(ii) - S(ji)i{u)) - (Sii2 - s li2 ) 
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xd»(i,k,(ji)) B»(l q ,(k(ji)) g ,(il) g ,2 q ) 4 2) (P(kUi)),P(ii)) 



+ 2 d° 3 (2, j, k) d§(5, l, (jk)) B° 4 (i q , i g , ((jk)l) g , %) 4 2 \ Pi ,p {ijk)l) ) 

-\ f 3 (i,j, k) 4(2, 1, (jk)) B° 4 (l q , (if) (Ujkj) %) Jf\p (h3) ,p um ) 



~ A° 3 (l,j, 2) 4(2, 1, k) B° 4 (l q , i g , (lk) g , %) 4 2 \ Ph p m ) 



(S 2 j(jk) - S2j(l(jk))) - (%)j(jfc) - S (ij)j((jk)l)) - {Slj2 - S lj2 ) 



x4(2,l,(jk)) Bl(\ q ,(ij) g ,((jk)l) g ,2 q ) 4 2 \p{i^P((jk)i)) 



+\ 4(2, I, k) 4(%j, (Ik)) B° 4 (l q , i g , ((lk)j) g , %) jf\p mj) ,Pi) 



-- f 3 (i, I, k) 4(2, j, (Ik)) B»(l q , (il) g , (j(lk)) g , 2 q ) 4 2 \p(ii),P(m))) 
-\ A° 3 (l, 1, 2) 4(2, j, k) B° 4 (l q ,i g , (jk) %) 4 2) (Pi,PUk)) 



{S2i(ik) ~ S 2 i(j(ik))) - {S(u)i((ik)j) - S(u)i(ik)) - (S112 - s ll2 ) 



x4(2,j,(lk)) B^(l q ,(il) g ,((lk)j) g ,2 q ) 4 2) (p (tl) ,p mj )) 



2 A° 3 (l, i, 2) A° 3 (l, k, 2) B° 4 (t q , I, \, %) 4 2) (p,, Pf ) 



+^ 4(1, 1, j) A° 3 (l, k, 2) Bl(l q , (ij) g , l g , %) J {2) (p & yPi) 

+\ 4(2, i, I) A° 3 (l, k, 2) B° 4 (l q ,j g , (U) g , %) 4 2) iPpPui 
1 

+ 2 



(S li2 - S m ) - (S li{ji) - - (S 2 i(u) - 



xA° 3 (l,k,2) B° 4 (l q ,(ji) g ,(il) g ,2 q ) 4 2 \p { J ir P ( u ) ) 



-\ A\(\, k, 2) A° 3 (l,i, 2) B2(f, J fl , 2 ? ) jf 

1. .11.: , . .n,: . 1. . : . , (2) 



+^ d°(i, fc, j) A° 3 (l, i, 2) lg, 2 q ) WiPfcyPi) 

+- d§(2, k, I) A° 3 (l,i, 2) (fcT) fl) 2 5 ) J 2 (2) (pj,P (fcl) ) 



+ - 



(#1*2 - -S'lfcs) - ( 5 lkCifc) - S Ik(jk)) ~ i S 2k(kl) - S 2 ~k(M)) 



xA° 3 (l,i,2) B° 4 (l q ,(jk) g ,(kl) g ,%) 4 2 \p i r k yP i ki ) ) 



-4(1, i,j) 4(2, 1, k) B° 4 (l q , (ij) g , (kl) g , %) 4 2 \p(i j ),P(kl)) 
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-4(ij,k)4(%iJ)B° 4 (t q ,(lk) g ,(ij) g ,%)4 2) (p m ,p m )j, (5.10) 

where Pc(i,j, k, I) denotes the set of cyclic permutations of the gluons. 

It is interesting to note the appearance of the sub-leading colour quark-antiquark 
antenna A\ and the accompanying A® antennae. The A\ antenna is introduced to remove 
spurious divergent behaviour of the D® antenna in the triple collinear limit. Specifically, 
the D® antenna contains the IR divergent limit involving colour-disconnected gluons, 

D°(i,j,k,l) Piji^,Q(x,y,z), (5.11) 

where the antenna tends to the QED-like triple collinear splitting function. This divergent 
limit has no analogue in the leading colour physical matrix elements and must be removed. 
In QCD the quark resides in the fundamental representation and so always terminates a 
string of gluons which lie to the right (left) of the quark (antiquark) in the colour ordering, 
e-g-, 

Mn(' " " ]lq,ig,jg,k g ,--- , l g , m g , n g , 2^; ■ ■ ■ ). (5.12) 

The D4 antenna is derived from the decay of a neutralino to a gluino and gluons, \ ~~ ^ 9999-, 
where the gluino plays the role of the quark. In contrast to the quarks, the gluino resides 
in the adjoint representation and so can sit anywhere in the colour ordering, including 
configurations with gluons either side of the gluino. For example, the D® is derived from the 
squared partial amplitude M®(i g ,j g , k g , l g ) which is symmetric under cyclic permutations of 
its partons. The cyclic symmetry of this matrix element gives rise to the spurious divergent 
limit described in Eq. (5.11). 

In any QCD calculation involving a quark-antiquark pair, the subtraction term will 
contain two D® antennae, one for the quark and the other for the antiquark, each contain- 
ing a spurious limit of the kind in Eq. (5.11). e.g., for the matrix element in Eq. (5.12), 
the subtraction term contains the antennae 3 , ^(l, k) and D °(2, n, m, I). The subtrac- 
tion term is constructed for a block of colour orderings related by cyclic permutations of 
gluons. Using this fact, the subtraction term will always contain the antennae 1)4(1, k) 
and 1)4(2, k, where the arguments of the second antenna are obtained by cyclically 
permuting the gluons in Eq. (5.12). Both of these antennae contain spurious limits of the 
kind in Eq. (5.11), involving gluons i and k and the antenna's fermion. The spurious limits 
of both D® antennae can therefore be removed using an A®(l,i, k, 2) antenna, as is done 
in Eq. (5.10). 

Using D® antennae to isolate divergences associated with the quark and antiquark 
endpoints of a quark string, forces us to construct the subtraction term for a block of colour 
orderings related by cyclic permutations of gluons. Exploiting the cyclic permutations, the 
spurious limits of the D® antennae can be systematically removed using an A® antenna. 
It is noted that this argument holds no matter how many gluons are present in the quark 

3 Here the quark-antiquark pair is in the initial state (IF D4 antennae) for consistency with the process 
considered in this section but the argument is generic to FF, IF and II antennae. 
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(a) 




(b) 
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(c) 



(d) 



Figure 6: Plots displaying the convergence of the subtraction term to the physical matrix element 
in various unresolved limits. The green data is furthest from the singular configuration with the 
blue data closer to the singular region and the red data the closest, (a) Final-state triple collinear 
limit between partons i,j,k such that x — Sijk/s\2- (b) Triple collinear limit between final-state 
partons i, j and initial-state parton 1, such that x = Suj/s^- (c) Double collinear limit between 
final-state partons i,j and the initial-final pair 1, k such that x = Sy/si2, y = s±k/si2- (d) Double 
soft limit for soft partons i,j such that x — (si2 — Sr/)/si2- 



string and the technique is not a special case of the six parton scattering process considered 
here. 

With the modifications due to the spurious limits of the quark-gluon antennae, the dou- 
ble real subtraction term fits the general structure derived in Sec. 3. Due to the non-trivial 
factorisation of the four-parton final-final antennae and the large sum of permutations, it 
is not straightforward to show analytically that this subtraction term correctly mimics the 
IR divergent behaviour of the physical cross section. However, to demonstrate its validity, 
the subtraction term has been tested numerically. 

For each IR divergent configuration, a set of momenta were generated using RAMB0 [70] 
such that the momenta fulfil a set of constraints defining the unresolved configuration. In 
this configuration, the ratio of the physical matrix element to the subtraction term we 



4 The number of gluons should be more than two. In the case of two gluons an A% antenna is used instead 
of the D*i and the issue does not arise. 
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Single collinear limit for qqbar -> gggg 
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Figure 7: Distributions of R without azimuthal angular rotations for single collinear limits: final- 
final collinear limit for (a) gg — > gggg and (b) qq — > gggg and initial-final gluon-gluon collinear limit 
for (c) gg — ¥ gggg and (d) qq — > gggg. The parameter x controls how closely the singular limit is 
approached, close (dark green), closer (blue) and closest (red). 



compute, 

T S 

NNLO 



R = "^ jVL ° . (5.13) 
do"? 



This procedure is then repeated for 10,000 different phase space points in the unresolved 
configuration defined by the constraints. The constraints are then tightened to force the 
phase space points closer to the unresolved singular point and the ratio calculated for 
another 10,000 points. The procedure is repeated once more for a set of points even closer 
to the singular point and the histogram of ratios for the three sets of constraints plotted. 

A selection of plots from four different unresolved configurations is shown in Fig. 6. 
In each plot we see that the distribution of R becomes more sharply peaked around one 
as the unresolved singular limit is approached, i.e., as the parameter x gets smaller. This 
provides statistical evidence for the convergence of the subtraction term to the physical 
cross section in the IR divergent limits. 

Testing the subtraction term in this way also gives a clear demonstration of the pres- 
ence of azimuthal correlations in the single collinear limits. The origin and solution to this 
problem have been discussed in [9,48,67], where it was shown that the azimuthal correla- 
tions which spoil the subtraction in the single collinear limit arise from splitting functions 
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where a parent gluon splits into two daughter gluons or a quark-antiquark pair (and the 
various crossings of this configuration to include initial-state partons). The qq — > gggg 
process contains the same final state as the gg — > gggg process, studied in [57]. There- 
fore by turning off the angular rotations when computing the ratio R, both subtraction 
terms should display similar behaviour when testing the single collinear limit between two 
final-state partons. This is seen in Fig. 7(a) and (b), where the subtraction terms for both 
processes display a broad peak, characteristic of the presence of uncompensated angular 
terms. An interesting comparison between the two processes is seen in the initial-final 
collinear limit. The collinear limit between an initial-state gluon and a final-state gluon 
does display azimuthal correlations while the collinear limit between an initial-state quark 
and a final-state gluon does not. We see that the limit involving the initial-state quark 
shown in Fig. 7(d) is sharply peaked while the corresponding limit in the purely gluonic 
process seen in Fig. 7(c) displays a broad peak because of the uncompensated azimuthal 
correlations. 



5.5 Real-virtual subtraction term, d&J^ NLO 

Following the discussion in Sec. 3, the real-virtual subtraction term is given by, 

1 dxi dx2 



dVqq,NNLO = ^5 J d$ 3 (P3, • ■ ■ , Pb\ ZlPl , X 2 p 2 ) ^ 



^(1<2' igiigi kg, 2q) B§(l q ,ig, jg, k g ,2 q ) ' (pi,Pj,Pk) 



X\ x 2 

(3), 



P(i,j,k) 



+ 4(i,i,j) 



+ J { l\l q ,{tj) q ,k g ,2 q ) Bl{l q ,{ij) q ,k g ,\) 



B\(l q , (ij) q , kg, 2 q )S(l - Xl )6(l - x 2 ) 



J { 2 2 \PipPk) 



+ f!(i,j,k) 



+ J { l\l q ,{i3) q ,{jk) q ,2 q ) Bl(l q ,(ij) q ,(jk) q ,2 q ) 



B\{l q , (ij) q , (jk) q , 2 q )5(l - Xl )S(l - x 2 ) 



j 2 \pi 3 ,p~ k ) 



+ 4(2, k,j) 



+ J^\l q ,ig,{kj) q ,2 q ) Bl{\ q ,i g ,{kj) q ,2 q ) 



Bl(l q ,ig,(kj) q ,%)6(l - x 1 )6(l - x 2 ) 



J?\pi,P~ k ) 



+ 



di(i,i,j)s(i - Xl )s(i - x 2 ) + j^\l g ,i g ,j g ) d° 3 (t,i, j 



-2 4(U,j) J ( 2 \l q ,(ij)) 



Bl(l q ,(ij) q ,k g ,2 q ) 4 2 \ Prj , Pk ) 



+ 



fl(i,j,k)5(l - Xl )S(l - x 2 ) + J { i\i g ,j a ,k g ) fl(i,j,k) 



-2fl(i,j,k) J { 2 l) ((ij) q ,(jk) 



Bt(l q ,(ij) g ,(jk) g ,2 q ) 4 2) (p rj ,P~ k 



+ 



4(2, k, j)5(l - Xl )5(l - x 2 ) + jg } (2„ k g ,j g ) 4i% k,j) 
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-2 4(2, k,j) J { 2 l) (2 q ,(kj) q ) 



B° 4 (l q ,i g ,(kj) g ,%) 4 2) (pi,p~) 



A£(l, i, 2)8(1 - Xl )5(l - x 2 ) + jf(l q , i g , 2 q ) A»(l, i, 2) 



-A§(l,i,2) J { 2 \l q ,% 



Bl(l q ,j g ,k g ,2 q ) 4 2) (p-, P - k ) 



+ & log (j^J rfg(l,i,j)(5(l - xi)*(l - x 2 ) B° 4 (I q , (ij) g ,k g ,2 q )J {2) (p Tj ,p k ) 



+ b log ( ^- ) /3 (i, j, fc)<5(l - si)<5(l - x 2 ) Bi(l q , (ij) (jk) g , 2 q )J {2 \p r p~ k ) 



+ fro log 



\l«2fcj| 



d§(2, fc, j)<5(1 - xx)5(l - x 2 ) Bl(l q , i g , (jk) g , %)jf\pi,p~) 



+ 



1 



\vl q (s m ) - \vl q ( Slj ) - \^(sr^ k ) + \^(s jk ) + Al qq (s V2 ) - Al qq (s r2 ) 



+5(1 - X1 )S(1 - x 2 ) (S(s^ )k , s jk , x^ )k jk ) - S(s jk , s jk , 1) - S( 8l jrr y s jk , x^ jk ) 
+S(si j ,s jk ,xi jdk )j 



1 

+ 2 



B 4 (I q ,(ij) g ,kg,2 q ) J ( 2 \p Tr P k ) 



-l^(s ik )+S(l- Xl )S(l-x 2 )( - 5(.s 



(ij)(jk)> Sik i x (ij)(jk),ik> 



3" • 5V (ij){jk)> 

+ S(s ik ,s ik ,l) 



+ ' S ( S ^ij) ,Sik, X l(ij),ik^ 5 ( S li> S iki x li,ik) + 5 ( S 2(jfc)' Sifc ' X 2(jk),ik) 



S( s 2k, s ik,X 2kyik )j 



fi(i,j,k) B° 4 (l q ,(ij) g ,(jk) g ,2 q ) J { 2 2) (p^p~ k ) 



+ 



Hh^)) " ^3° " ^3°(-^) + ^(%) + Ai qq (8 l2 ) - Al qq (.S- l2) 



+ 



3" 6 ' 3 

+ 5(1 - Xl)S(l - X 2 )(S(s.^ ky Si j ,X.^ k) ..) - S(8ij,8ij, 1) 
~^( S 2(&y Sij,X 2{kj),i^ + 5 ( S 2j) s ij, x 2j,ij)) 

x d§(2, B° 4 (l q ,i g ,(kj) g ,%) J {2) ( Pi ,p~) 



-2 P 3, g ( s 2fc) + S(l - xi)<5(l - a:2)^(s I 2,s^,x I 2^) - S(s- 12 ,s jk , X - l2Jk ) 
S(si~-,s-- k ,xi~-~- k ) + S(sij, Sj k , xijjk) - S( s %, s ]h x 2k~j 



+S(s 2k , s jk,X 2k j k ) S j 



jki Jj 2k,jk) 

Al(l,i,2) B 4 (l q ,jg,kg,%) A^iPj^k) 



(5.14) 



The integrated antenna strings used in this subtraction term are given by, 

■^5 ^(^?> igijgi kg, ^g) = \^-g> ig) ^2 \igijg) ^2 \jg> kg) + J 2 \^ q , kg), 



JbO-qi igiJg) — J 2 \^qi ig) + J\ '{igiJg) + J 2 '(2 q ,j g ), 



r(l), 
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39i kg) ~ J 2 \igi3g) + J 2 \jg> k g ) + J 2 \ig, k g ), 
jl D(2q,igJg) = jl D0-qiigi3g)i (1 +* 2, zi« z 2 ), 

Ji(i 9 ,i fl ,2,-) = 4 1) (I 9 ,2,-)- (5-15) 

After integration only those terms introduced at the real-virtual level are passed down to 
the double virtual subtraction term. 

5.6 Double virtual subtraction term, ^Ytnlo 

Following the discussion in Sec. 3.3, the double virtual subtraction term is constructed 
from the remaining terms inherited from the double real and the real-virtual subtraction 
terms, in addition to the double virtual mass factorisation contribution. The resulting 
double virtual subtraction term can be written in terms of single and double unresolved 
integrated antenna strings, 



MqqNNLO = -JV 4 a$ 2 (P3,P4;Pl,P2) / — XT 2_> 



J^O-qiigi jgi^q) (^4 0-qi 3gi 2g) ~ B^(l q ,ig, j g ,2 q 
+ 2 \^-qi igi3g-> 2q) ® J4 \^-q, ig, 3gi 2q) B2(l q ,ig, jg,2q) 

+ jf (i^,.^-) B2(i„»p,i fl ,2,-)} 4 2) (puPj)- (5-i6) 

The single unresolved integrated antenna string is defined in Eq. (4.12). The new ingredient 

(2) 

at the double virtual level is the double unresolved integrated antenna string, J\ , which 
displays the structure discussed in Sec. 3.3. In Eq. (4.12) the single unresolved integrated 

antenna string was written as a sum over integrated dipoles. The analogous dipole-like 

(2) 

formula for J \ is given by, 

jf\t q ,i g ,j g ,2 q ) = J { 2 2 \l g ,i g )+ J { 2 2 \ig,j 9 )+ J { 2 2) {3gX)- J?\hM (5-17) 
where the two-parton double unresolved integrated antenna strings are given by: 

J?\i q ,i g ) = \vl q ( sil ) + \vl q (s- u ) + e vl q (s- u ) 

- \ Kq(s- U ) ® Vl q { Sll )] ( Zl ) - T q l\ Zl )6(l - Z2), (5.18) 

J?\ig,3g) = \n{s i3 ) + \H{s l3 ) + | (^) " V 3 °(%) 

" (5-19) 
jf(i 9 ,2^-) = jf(2,,i 9 ). (5.20) 

An integrated dipole, J 2 \i, j) is associated with a power of (|sjj|)~^ e , thereby matching 
specific terms in the Catani representation of the one- and two-loop matrix elements. 
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However, for the leading colour contribution to qq — > gg jets, we do not expect a singular 
contribution that depends on s\2- 

However, the D\ antenna contain infrared singular limits that are not present in the 
matrix elements. These singularities are compensated by the A^ qq contribution in dafj NLO 
and by the A\ qq and and A® qq - x A® terms in daJj NL0 . Terms of the form (V® q x A®) and 
(.A3 qq x d®) which were introduced at the real-virtual level, cancel after integration. The 
integrated forms of these spurious contributions reappear in da^ NLO and are collected in 



7 "(2) 
J 2 1 



J?\hA) = - \[A,d s i2) ® Al qq (s u )](z 1 ,z 2 ), (5.21) 

which is proportional to (|si2l)~ 2e - At first sight, the presence of J 2 (l q ,2 q ) appears to 
be in conflict with our earlier statement about the absence of singularities proportional to 
si2- However, the leading singularity in J 2 is 0(l/e), 

Voles (eJ 2 2 \l q ,2 q )) = 0, (5.22) 

so that, after expansion of the ( | ^12 1 ) 2e factor, the singularity does not depend on sjg 
while the remaining finite terms can and do depend on ln(sj2"). 

5.7 Comparison with the IR structure of the e + e~ — > 3 jet leading colour double 
virtual subtraction term 

It is interesting to re-examine the form of the double virtual subtraction term for the 
e + e~ —> 3j process [48]. Although the e + e~ annihilation calculation involves only final- 
state partons, the integrated form of the subtraction terms is expected to closely follow the 
qq —7- gg double virtual subtraction term up to mass factorisation contributions and crossed 
antenna. The double virtual subtraction term as presented in [48] can be reformulated in 
terms of final-state integrated antenna strings to fit the double virtual structure discussed 
in Sec. 3, so that, for example, at leading order in the number of colours, 

da% 2 = A3 N 2 d<S> 3 ( Pl ,p 2 ,P3,q 2 ) 

4\U 3 9 , 2 q ) UfUl q , 3 g , 2 q ) - b f M 3 °(l ? , 3 g , 2 q ] 

+ l - j£\l q ,3 g ,2 q )®jW(l q ,3 g ,2 q ) Ml(l q ,3 g ,2 q ) 

+ 4 2) (l„3 fl ,2 ? ) M 3 (l ? ,3 ff ,2 ? -)} 4 3) (pi,p 2 ,P3), (5.23) 

where M£(l„ 3 g , 2 q ) and M 3 1 (l g , 3 g , 2 q ) are the tree-level and one-loop squared partial 
amplitudes for the process e + e _ — > 7* /Z — > qgq and q 2 is the virtuality of the colourless 
boson. The factor A3 carries all overall colour factors, QCD coupling and non-QCD factors, 

A3 = 47ra ^e 2 q g 2 {N 2 - I). (5.24) 
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The final-state single unresolved integrated antenna string is given by, 

^3 \lqi ^g, 2g) = J 2 \lq, 3 fl ) + J2 \3g, 2g), (5.25) 
4 1) (3 s ,2 9 -) = j( 1) (2„3 s ) (5.26) 

The double unresolved integrated antenna string for the relevant final-state partons is given 
by, 

jf\l q ,3 g ,2 q ) = 4 2) (l g ,3 9 ) + 4 2) (3 9 ,2 g -)- J?\l q ,2,), (5.27) 

where 

J?\l q ,3 g ) = l -Vl{ Sl ,) + l -V\( Sl ,) 

+ h f e {^) ^3(^13) " ^(«is) ® ©§(-13), (5-28) 

4 2) (3 fl ,2 9 -) = Jf(2„3 s ), (5.29) 

In this example, all integrated antennae are of final- final type and so carry all z\ and z 2 
dependence though (^-functions, rendering all convolutions and integration over z\ and Z2 
trivial. As expected, the relevant integrated antenna strings are composed of two-particle 

integrated dipoles, and J^. 

(2) 

Once again, J 2 is present to remove the spurious singularities generated by the D 
antennae, 

J?\l q ,2 q ) = \a%s 12 ) + A\{s 12 )-\aI{s 12 )®AI(s 12 ). (5.30) 
It has exactly the same origin, structure and singular behaviour as Eq. (5.21). 

6. Conclusions and discussion 

In this paper, we made a detailed investigation of the infrared structure of NLO and NNLO 
pertubative calculations where the antenna subtraction method is used to isolate the in- 
frared singularities present in the radiative contributions to the cross section. At NLO, 
there is a clear and simple relationship between the unintegrated antenna used to subtract 
the implicit singularity present in the real radiation contribution and the (singular) inte- 
grated antenna that cancels the explicit poles in the virtual contribution. At NNLO, the 
subtraction takes place across the double real, real-virtual and double virtual contributions 
and is significantly more complicated than at NLO. However, the implicit divergence of 
the double real and real-virtual contributions is captured by distinct blocks of subtraction 
terms which follow a predictive structure based on the colour ordering of the matrix ele- 
ments under consideration. Fig. 3 shows explicitly how the integrated and unintegrated 
subtraction terms conspire to cancel the explcit divergences in one- and two-loop matrix 
elements. 

Each unintegrated antenna involves two hard radiators that, after integration over the 
unresolved phase space, appear as a pair of colour connected final state particles, I and 
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J, in the virtual subtraction term. The singularities associated with a integrated antenna 
are proportional to (\su\)~ ie and this gives a natural linkage to the one- and two- loop IR 
singularity operators and -T^(e) proposed by Catani [64]. We denote the collection 

of integrated antennae for a particular pair of colour connected particles (plus the relevant 
mass factorisation contribution) as an integrated dipole. At one-loop, is related to an 

integrated three-particle tree-level antenna and describes the unresolved radiation between 

(2) 

two colour connected particles. Similarly, at two-loops, J \ fulfills the same role for double 
unresolved radiation and involves integrated four-particle tree-level antennae, three-particle 
one- loop antennae and products of three-parton tree- level antennae. Mass factorisation 
contributions are absorbed into the integrated dipoles so that the deepest poles contribute 
only in the soft region. 

The IR structure of the two-loop contribution is determined by integrated antenna 
strings which are simply formed from J 2 ^ and ■ For a given colour ordering of a 
particular n-particle process, Jn ^ and J n ^ is simply a sum over integrated dipoles that 
involve colour connected particles. Equivalently, the pole structure is identified as a sum 
of dipole-like terms proportional to (\sij\)~ et (for i = 1,2) that link the colour connected 
particles i and j. There is one caveat: When the integrated dipole involves the and 
D\ antenna, there is a well defined correction factor J 2 that removes the spurious colour 
connections between the quark and antiquark. 

The integrated antenna strings contribute to the real-virtual and double virtual sub- 
traction terms in a very precise manner, again as illustrated in Fig. 3, leading to the 
NNLO master formula of Eq. (3.62) for the double virtual subtraction term. In this form, 
there is a direct correspondence with Catani's one- and two-loop factorisation formulae 
and the explicit pole cancellation against the double virtual contribution is particularly 
straightforward. 

The form of J$ and jffl in terms of integrated antennae imposes a particular structure 
on the unintegrated antennae that make up the real radiation subtraction terms. Under- 
standing the explicit pole structure of virtual amplitudes in terms of integrated antenna 
strings gives a direct connection between the block structure of the unintegrated subtrac- 
tion terms and the explicit pole structure of virtual contributions, thereby simplifying the 
construction of the double real, real-virtual and double virtual subtraction terms. 

We anticipate that the structuring of the NNLO subtraction terms da^ NLO , daJj NLO 
and dcijy^Q, presented here will be of great help towards the evaluation of the NNLO cor- 
rections to a number of important processes involving jets at hadron-hadron and electron- 
hadron colliders such as pp — > 2 jets, pp —> H+ jet, pp — > V+ jet and ep — > 2(+l) jets 
where the experimental interpretation is limited by the theoretical uncertainties. In each 
of these cases, the ultimate goal is a full parton-level Monte Carlo implementation. 
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A. Collinear splitting kernels 
A.l Leading order splitting kernels 

The splitting kernels present in the definition of the mass factorisation terms (2.25), (3.2) 
and (3.3) contain colour factors which can distribute the colour ordered splitting kernels 
across multiple orders in the colour decomposition. 

The leading order splitting kernels can be decomposed into colour factors N and N F . 
Expanding the splitting kernels into this set of colour factors defines a set of colour stripped 
splitting kernels: 



N 2 



1 



N 



Pqg{ x )i 



Pqq{x), 
Pgq(x), 



P° gg (x) = N p° gg (x) + N F p° mF (x). 
In this decomposition, the colour ordered splitting kernels are given by [43] , 



(A.l) 



Pqq( X 



P gq(x 



Pqg( X 



p° gg (x 
P°ggA x 



^ . . (1 + x) 3 . 
V (x)- K -^ + -5(l-x), 



x 

1 



1 + 



X + X 



2' 

2 



(A.2) 



2V (x) + - - 4 + 2x - 2x 2 + b 5(l - x), 

x 

bo,FS(l - x), 

where the distributions V n {x) are defined in terms of plus-distributions, 

m , , /ln"(l-x)\ lk , 

The mass factorisation kernels are classified using a similar notation to the splitting kernels, 

' N 2 — V 



r«(x) 



N 



r£ } (x) = 



N 



1 



r (i)( 
x qg v 



N 

~^qg\ X )i 



Tqq( X )> 

r«(x), 



rg(x) 



NTW(x) + N f tW (x) 



99- 



(A.4) 
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where the mass factorisation kernels in the MS scheme are related to the splitting kernels 

by [58], 

rg>(s) = --P%(x). (A.5) 
A. 2 Next-to-leading order splitting kernels 

The one-loop quark-quark and quark-antiquark splitting kernels, P qiQj and P qiqj , contain 

s 1 

a non-trivial flavour structure and are classified according to flavour singlet, P ^ , and 
flavour non-singlet, P^j 1 contributions [71]. 

^ > '?i'?j( X ) = ^V^qq ^qq ( X ^> 

P\ iq .(x) = SijP^ix) + P%(x), (A.6) 

where P| : (x) = P S qq \x). 

For the purposes of the DGLAP evolution, the parton distributions are decomposed 
into a flavour basis of non-singlet, f n (x,/j, F ), and singlet, f s (x,/jL F ), combinations, where 
x denotes the fraction of the hadron momentum carried by the parton and \j? f is the 
factorisation scale. Suppressing the x and n F dependence of the distributions, the flavour 
asymmetry distributions are given by, 

fn,ij = (fi±fi)-(fj±fj)- (A.7) 
The total valence distribution is given by, 

Np 

f V n=T,(fi-fi), (A.8) 
i=\ 

and the singlet distribution, which mixes with the gluon distribution during DGLAP evo- 
lution, is given by, 

Np 

fs = J2(fi + fi). (A.9) 
i=i 

Each distribution is evolved between momentum scales using the corresponding DGLAP 
evolution kernel. The flavour asymmetry, total valence and sea evolution kernels are de- 
noted Pn(x), Pn( x ) an d P s (x) respectively. These kernels are related to the singlet and 
non-singlet splitting kernels in the following way: 

Pt(x) = pV(x)±pV(x), 
Pl(x) = P~(x) + P s n (x), 

P s (x) = P+(x) + P s s (x), (A.10) 

where the non-singlet kernel, given by difference between the total valence P v n {x) and the 
"minus" -P^(x) evolution kernels is, 

P s n (x) = N F (P s qq (x) - P s qq (x)), (A.ll) 
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and the "pure singlet" kernel, given by the difference of the sea and "plus" evolution kernels 
is, 



P s s (x) = N F (P s qq (x) + P s qq (x)). 
The singlet and non-singlet splitting kernels have the perturbative expansions: 



(A.12) 



P v qq {x) = P° q Jx) + 



2^ 



P v qq \x) + 0{ai), 



qq 



2tt ) 
2vr / 



)P^(x) + 0(a 2 ), 
)P s qq 1 (x) + 0(a 2 ), 



)P S qq 1 (x) + 0(a 2 ). 



(A.13) 



At one-loop, the relation P q ^(x) = P qq (x) simplifies the perturbative expansion of 
Eq. (A.10): 



Pt l (x) 
P^\x) 



pV q \x)±P^(x), 
P+'\x) + Pf(x), 



V,l, 



where the one-loop pure singlet kernel is given by, 

Pf(x) = 2N F P s qq \x). 



(A. 14) 



(A.15) 



The one- loop non-singlet splitting kernels, P^. q {x) and P q ^(x), are presented explicitly in 



Refs. [71,72], and the singlet evolution kernel, P],(x), in Ref. [71]. Rearranging for P q ' q (x) 
in Eq.(A.14) yields the quark-quark singlet splitting kernel, 



P S qq 1 (x) = ^-{P\(x)-Pt: 1 (x)). 



(A.16) 



This function can be derived from the forms of P^i 1 (a;), P ' q (x) and P\{x) in Ref. [71] 



using Eq. (A.16) or found explicitly in Ref. [73]. Given the explicit form of -P^ 1 (x), P^(x) 
and Pq^x), the full one-loop splitting kernels, P qq (x) and P qq (x) are easily obtained. 
The one-loop splitting kernels can be decomposed into the N, Np colour factors: 



P\ q {*) 
P gg -(x) 
P\ Q (x) 
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N 2 - 
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N P\q( X ) + Pqq( X ) + l^T + N F PqqA X ) 



■qq 



qq, 



P\q{ X ) + Pqq{x) 



N 



p\q{x), 
p\q( x )i 
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N 2 -l 
N 



N P] q ( X ) + PU X ) + N F PqqA X ) 



N 1 9q 
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P\ g {x) = Np\ g (x) + ^p\ g {x 



P] g ( x ) = N 2 p gg {x) + NN F p] g ^ F (x) + ^ p] g , F { x ) 



(A.17) 



At NNLO, two new ingredients are present in the mass factorisation terms: convolutions 
of one-loop mass factorisation kernels, either with other mass factorisation kernels or inte- 
grated antenna functions: 

[ r Sf ® rffe] *a) = J dxxd yi rg ) (x 1 )r«(y 1 )«5(^ 1 - x iyi )S(i - z 2 ), 

[rg } ® *3 ] (a; z 1; z 2 ) = J dzidj/i rg^Xi)^ ^; y 1 ,z 2 )S(z 1 - x m ), (A.18) 



and the reduced two-loop mass factorisation kernel, defined as, 



=1(2)- 



1 



2c 



u 1 



(A.19) 



Decomposing the reduced two-loop mass factorisation kernels into the N N F colour factors 
yields: 
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=;(2) 



If (2) ■ 
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T^(2) 



N 2 T g 2 g \x) + A^rf^x) + ^f^ F (x) + N F T g 2 g \ F2 (x) 



(A.20) 



where the explicit form of each colour stripped function can be simply inferred from the 
colour decomposition of the leading order and one-loop splitting functions in Eqs. (A. 2) 
and (A.17), together with the definition of the reduced two-loop kernel in Eq. (A.19), to 
be: 



r q 2 q \x) 
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B. Integrated large angle soft term: Initial-final configuration 

In this appendix, we give an expression for the integrated large angle soft terms when the 
primary mapping is of the initial-final form. This is an integral of the soft antenna function, 

2 Sac 



S, 



ajc 
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over the unresolved initial-final phase space for the mapping (i,j,k) — > (I,K) defined 
by [44], 



S IF (Sac, SIR, Uac,iK) 
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(B.l) 



After integration over the antenna phase space [44] and expanding in distributions," 
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